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INTRODUCTION 


1 . 


The  advent  of  very  large  transport  aircraft  with  the  powerful 
wake  vortices  which  they  produce  has  introduced  a serious  threat  to 
smaller  aircraft . ^ ^ ^ ^ ^ Numerous  upsets  and  disturbances,  because 
of  vortex  encounters,  have  been  documented,  and  some  serious  aircraft 
accidents  have  been  traced  to  this  source.  The  problem  is  most  acute 
on  landing  approach  and  during  climbout  after  takeoff  when  flying 
speeds  are  low  and  precise  path  control  is  important. 

An  extensive  theoretical  and  experimental  vortex  research  program 
is  being  pursued  by  the  U.S.  Department  of  Transportation,  Transporta- 
tion Systems  Center  (DOT/TSC)  . ^ Results  of  this  program  have  had  a 
major  impact  on  the  basic  understanding  of  trailing  vortices.  At  a 
number  of  large  airports,  TSC  has  demonstrated  an  effective  method  of 
detecting  vortices  and  locating  their  lateral  positions  near  runways. 
The  method  uses  ground-wind  sensors  in  the  form  of  propeller  anemo- 
meters, placed  in  lines  perpendicular  to  the  approach  path  of  aircraft. 
Software  developed  by  TSC  to  postprocess  data  from  these  anemometers 
produces  a graphical  representation  of  probable  vortex  positions. 
However,  this  software  is  only  partially  automated  and  requires  sub- 
stantial manual  pattern  recognition  and  manual  processing  to  reduce 
the  data  completely. 

This  report  documents  efforts  at  The  Charles  Stark  Draper 
Laboratory,  Inc.,  under  contract  to  TSC,  to  automate  fully  the  data- 
reduction  task.  The  software  system  which  has  been  developed  processes 
propeller-anemometer  data,  and  automatically  tracks  and  records  vortex 
position  as  a function  of  time.  In  addition  to  vortex  tracking,  the 
system  also  provides  automatic  sensor-failure  detection  and  identifi- 
cation by  using  the  multiple  redundancy  of  sensors  in  the  sensor  lines. 
With  minimal  modification,  the  software  can  be  readily  adapted  to 
provide  real-time  tracking  in  planned  future  wake-vortex  avoidance 
systems . 
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In  the  following  sections,  the  relevant  aspects  of  the  software- 
development  effort  are  described.  A detailed  description  of  the  basic 
vortex-tracking  algorithm  is  presented.  Accompanying  this  are  des- 
scriptions  of  the  techniques  used  to  detect  the  onset  of  useful 
tracking  data,  to  determine  when  useful  data  are  no  longer  available, 
and  to  eliminate  siurious  bad  data  when  they  occur  during  an  otherwise 
useful  tracking  period.  Results  of  testing  the  algorithms,  using 
actual  sensor  data  gathered  by  the  TSC  system  installed  at  London's 

g 

Heathrow  International  Airport,  are  explained  in  detail,  and  the 
character  of  the  tracking  algorithms  is  discussed.  The  methods 
developed  for  sensor-failure  detection  and  identification  are  presented, 
and  the  results  obtained  using  the  Heathrow  data  are  examined.  Computer 
requirements  to  implement  the  software  in  a fielded  system  are 
summarized.  Finally,  a tentative  plan  for  the  future  development  of 
a simulation  of  a fielded  vortex-tracking  system  is  given. 
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2.  VORTEX-TRACKER  DEVELOPMENT 


The  vortex-tracker  development  is  based  on  optimal  estimation 

9 

theory,  using  appropriate  mathematical  models  to  describe  vortex 
behavior.  These  models  are  the  simplest  possible  descriptions  which 
are  consistent  with  the  tracking  requirements  and  embody  only  those 
elements  which  are  essential  to  attain  desired  tracking  performance. 

The  following  sections  describe  the  models  and  the  vortex  tracking 
algorithms  which  were  derived  from  them. 

2 . 1 VORTEX  DYNAMICS 

At  distances  of  five  or  more  wingspans  behind  an  aircraft,  the 
flow  field  generated  by  the  aircraft  is  approximated  by  two  line  vortices, 
one  trailing  from  each  wing  tip.  These  vortices  are  of  equal  and 
opposite  strength  and  are  separated  by  a distance  which  is  approximately 
equal  to  the  aircraft  wingspan.  The  flow  field  generated  by  each  vortex 
is  a rotational  motion  in  which  the  fluid  flows  in  concentric  circles 
about  the  vortex  center,  so  that  the  velocity  is  always  tangential.  At 
large  distances  from  the  center,  the  flow  velocity  is  proportional  to 
the  vortex  circulation  and  inversely  proportional  to  the  distance  from 
the  center.  The  usual  expression  for  tangential  velocity  is 


where 


Vg  = tangential  velocity, 
r i vortex  circulation,  and 
r : radial  distance  from  vortex  center. 
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From  Eq.  (2,1),  it  is  clear  that  as  the  center  of  the  vortex  is 
approached,  the  tangential  velocity  increases  without  bound.  Laboratory 
and  full-scale  measurements  have  shown  that  in  an  actual  vortex  there 
is  a region  near  the  center  where  the  i model  is  invalid  because 
the  flow  velocity  tends  to  zero  as  the  center  is  approached.  This 
region  is  defined  as  the  vortex  core  and  for  most  aircraft  of  interest, 
the  core  of  a newly  formed  vortex  is  smaller  than  about  10  feet  in 
radius.*  There  is  limited  Icnowledge  of  the  details  of  flow  within  the 
core  and  a number  of  models  have  been  proposed.  These  models  generally 
provide  a tangential  velocity  field  with  velocity  increasing  monoton- 
ically  with  distance  from  the  center.  Parameters  are  chosen  so  that  the 
velocity  of  the  core  model  matches  the  velocity  of  the  i model  at  the 
radius  of  the  core  boundary.  The  largest  uncertainties  with  this  type 
of  model  are  usually  associated  with  judging  the  size  of  the  core. 

The  vortex  center  tends  to  move  in  space  according  to  the  velocity 
field  in  which  it  is  placed.  The  two  vortices  generated  by  an 
aircraft  interact  and  produce  vortex  translational  velocities  which 
are  orthogonal  to  the  line  between  their  centers.  Distances  between 
them  are  almost  always  many  times  the  core  size.  Hence,  the 
velocity  induced  by  each  vortex  on  the  other  is  determined  by  Eq.  (2.1) 
as 


where  b is  the  distance  between  the  two  vortices.  Thus,  if  the  two 
vortices  generated  by  an  aircraft  are  at  the  same  altitude  (i.e.,  the 
line  between  them  is  horizontal) , then  they  translate  downward  at 
the  velocity  given  in  Eq.  (2.2), 

In  a similar  manner,  if  two  vortices  of  equal  and  opposite 
strength  are  situated  one  above  the  other  and  separated  by  a distance 
d,  they  will  move  horizontally  at  a speed  r/2Tid.  Furthermore,  a 
horizontal  plane  of  symmetry  will  exist  half-way  between  the  two 
vortices  and  there  will  be  zero  induced  velocity  across  the  plane. 

As  the  lower  vortex  is  the  image  of  the  upper  vortex,  the  flow 
above  the  plane  is  seen  to  be  the  same  as  that  of  a single  vortex 


*A8  the  vortex  ages,  the  core  radius  may  grow  to  30  feet  or  more. 
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of  strength  T,  and  at  height  d/2  above  the  ground.  Hence,  the  ground 
effect  produces  a horizontal  velocity  on  the  vortex  according  to 


where  h is  the  altitude  of  the  vortex  (i.e.,  half  the  distance  to 


the  image) . 

If  there  is  wind,  then  the  vortices  also  move  according  to 
the  velocity  of  the  wind  field.  A steady  ambient  wind  produces  a 
steady  translation  of  the  vortices  at  the  ambient  wind  velocity. 

The  total  effect  of  the  interaction  of  the  vortices  with  each 
other,  influence  of  the  ground, and  influence  of  the  wind  is  additive. 
Interaction  of  the  two  vortices  with  each  other  produces  a downward 
motion.  When  the  vortices  pass  below  approximately  100  to  150  feet  of 
altitude,  depending  upon  vortex  strength,  the  ground  begins  to  affect 
their  motion  and  they  migrate  away  from  each  other.  A steady  wind 
tends  to  produce  motion  in  the  direction  of  the  wind.  The  influence 
of  the  three  effects  on  lateral  vortex  motions  is  illustrated  in 
Figure  2.1. 


PORT- 

PLANE 

-V  STARBOARD- 

VORTEX 

IMAGE 

• )-r  VORTEX 

V'  IMAGE 

FIGURE  2.1.  TYPICAL  VORTEX  PATHS  IN  THE 
PRESENCE  OF  STEADY  WIND 
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As  time  passes,  the  vortices  gradually  lose  their  strength  as 
the  result  of  a number  of  fluid  dynamic  mechanisms . ^ ^ Air 

turbulence  is  the  most  important  mechanism  promoting  vortex 
dissipation.  Small  scale  turbulence,  of  the  order  of  the  vortex 
core  size,  produces  a diffusion  of  the  vorticity  resulting  in  a 
sudden  "burst”  of  the  core  to  an  order  of  magnitude  larger  in  size. 

Rapid  decay  of  the  vortex  often  follows  the  burst  phenomenon. 

Turbulent  eddies  of  a few  hundred  feet  in  scale  length  tend  to  promote 
a phenomenon  studied  by  Crow  , which  distorts  the  line  vortices  into 
a series  of  ring  vortices.  Commonly, these  ring  vortices  deteriorate 
through  the  "burst"  mechanism  described  above.  Under  certain 
conditions,  however,  the  ring  vortices  may  persist  and  pose  a threat 
to  following  aircraft.  Very  large  scale  turbulence  tends  to  produce 
a general  motion  of  the  vortices  as  with  ambient  wind. 

2.2  SIMPLIFIED  VORTEX-TRANSPORT  MODEL 

The  information  for  use  in  tracking  the  vortices  consists  of 
the  outputs  of  ground-wind  anemometers.  Typical ly» they  are  placed 
a few  thousand  feet  before  the  aircraft  touchdown  point  and  in  lines 
perpendicular  to  the  runway  centerline.  Following  passage  of  an 
aircraft  over  the  sensor  lines,  the  trailing  vortices  descend  until 
they  begin  to  experience  the  ground  effect.  The  onset  of  ground  effect 
usually  coincides  with  the  time  at  which  the  vortices  begin  to  have 
a measurable  effect  on  the  sensors;  the  time  when  useful  sensor 
data  begin  to  become  available.  Once  in  ground  effect,  the  vortices 
separate  rather  rapidly,  and  their  interactions  effectively 
vanish.  They  tend  to  move  at  relatively  constant  velocity 
primarily  under  the  influence  of  ground  effect  and  wind.  Thus,  during 
the  period  when  useful  ground-wind  anemometer  measurements  are 
available,  the  downwind  vortex  travels  at  a nearly  constant  velocity 
which  is  the  sum  of  local  wind  and  ground  effect.  Simultaneously, 
the  upwind  vortex  also  moves  at  nearly  constant  velocity,  but  the 
wind  and  ground  effect  are  in  opposition,  and  hence  these  effects 
subtract  from  each  other  (see  Fig.  2.1). 

The  physical  character  of  the  vortex  tracking  problem,  as 
described  above,  lends  Itself  to  a simple  mathematical  model  for 
vortex  motion.  The  basic  assumption  is  that,  during  the  period  when  the 
vortex  produces  useful  sensor  information,  its  transport  velocity  is 
relatively  constant  and  is  dominated  by  ambient  wind  and  ground  effect. 
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Defininq  x as  the  horiEontal  distance  measured  alonq  the  line  of 
sensors  from  runway  centerline  to  the  vortex,  the  equation  of  vortex 
motion  is 

X = u + V,  (2.4) 

ambient  lateral  wind  velocity  alonq  the  sensor  line,  and 
the  sum  of  vortex  velocity  caused  by  ground  effect, 
the  difference  in  wind  velocity  from  ground  to 
vortex,  wind  shear,  and  other  unmodeled  effects. 

Typically ^ there  are  many  ground-wind  sensors  (i.e.,  of  the  order 
of  20) , so  that  the  ambient  qround-wind  velocity  u can  be  estimated  to  high 
accuracy.  Hence,  u is  assumed  to  be  a known  input  to  the  system 
because  errors  in  determining  u are  quite  small.  The  other  velocity 
variable  v is  intended  to  model  the  remaining  effects,  including  the 
difference  in  wind  from  ground  level  to  the  vortex,  ground  effect, 
and  a number  of  other  unknown,  unmodeled  or  random  effects.  The 
equation  for  v is 


where 


u 

V 


V = w , 


(2.5) 


y^j^ere  = white  process  noise. 

In  other  words, the  variable  v is  a Brownian  motion, and,  therefore, 
it  is  characteristically  an  aimless,  wandering  type  of  random  process. 
As  such,  it  is  used  to  account  for  a number  of  effects  wliich  are  quite 
random  in  nature,  or  for  which  there  is  very  limited  knowledge  in 
terms  of  accurate  modeling. 

2.3  VORTEX-INDUCED  VELOCITY  MODEL 

During  the  period  when  useful  data  are  available  from  the  ground- 
wind  sensors,  the  vortex  centers  are  many  core  radii  from  the  sensors. 
Hence,  the  i flow  model  is  appropriate  to  determine  the  component  of 
horizontal  velocity  measured  by  the  sensors.  Since  the  sensors  are  at 
approximately  ground  level,  ground  effect  is  important  and  must  be 
accounted  for.  Hence,  the  flow  field  is  induced  by  both  the  actual 
vortex  and  its  image.  The  situation  is  depicted  in  Figure  2.2.  A 
vortex  and  its  image  each  produce  circular  flow  fields  with  tangential 
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FIGURE  2.2  VORTEX  AND  GROUND-WIND  SENSOR  GEOMETRY 

velocity  es  given  in  Eg.  (2.1)  . At  the  ground  plcme.  the  vertical  coiT\ponents 
of  velocity  cancel  and  the  horizontal  components  add.  Thus,  if  there 
is  a ground-wind  sensor  at  the  point  d,  it  measures  a horizontal 
component  of  velocity  given  as 


v„  “ — sinQ  = 
m Tir 


hr 


nr 


Th 

ii(h^  + (x-d)^l 


(2.6) 


where  9 is  the  elevation  of  the  vortex  relative  to  the  ground,  as 
measured  from  the  sensor  position,  x is  the  lateral  position  of  the 
vortex^ d is  the  sensor  location, and  h is  vortex  altitude. 


From  Eg.  (2.6),  it  is  clear  that  the  effects  of  the  vortex  on 
the  sensor  measurements  are  determined  by  three  variables  t circulation. 


altitiide  , an>l  horizontal  distance  from  sensor  to  vortex.  Hence,  in 
theory,  simultaneous  measurements  from  three  anemometers  can 
determine  vortex  location  and  strength.  In  fact,  attempting  such  a 
simultaneous  solution  produces  rather  poor  estimates  of  vortex  loca- 
tion. Sensor  errors,  modeling  errors,  and  local  wind  all  tend  to 
degrade  the  measurements  to  the  point  where  it  is  necessary  to  use 
another  approach. 

2.4  GROUND-WIND  ANEMOMETER  PREPROCESSING 

The  ground-wind  anemometers,  which  provide  the  primary  tracking 
information,  are  grouped  in  so-called  sensor  lines.  Each  line  consists 
of  about  20  sensors  placed  on  a line  orthogonal  to  the  runway  center- 
line.  The  sensors  are  spaced  at  50-foot  intervals,  with  the  middle 
sensor  located  on  the  runway  centerline  and  typically  at  a distance 
of  from  1000  to  300  feet  before  the  runway  threshold.  Each 
sensor  measures  the  component  of  local  air  velocity  which  is  horizontal  and 
orthogonal  to  the  runway  centerline.  Positive  velocity  is  defined  as 
left  to  right  across  the  flight  direction. 

There  are  three  primary  elements  which  affect  the  sensor  output: 
vortex-induced  velocity,  local  wind, and  sensor-instrument  errors  or 
noise.  The  vortex-induced  velocity  component  is  discussed  in  Section 
2.3  and  is  given  in  Eq.  (2.6).  Local  wind  is  an  extremely 
variable  phenomenon,  depending  primarily  upon  meteorological  factors. 

For  the  purpose  of  this  discussion,  the  wind  effects  will  be  separated 
into  two  categories  depending  upon  the  scale  of  turbulence  in  the 
air.  One  category  consists  of  all  components  of  local  wind  which  are 
attributable  to  air  masses  of  a scale  which  is  large  compared  to  the 
length  of  the  sensor  line  (i.e.,  of  scale  length  greater  than  about 
1000  ft) . This  type  of  wind  involves  coherent  motion  of  large  masses 
of  air, and,  hence,  it  tends  to  affect  all  sensors  in  a sensor  line  in 
the  same  way.  A steady  ambient  wind  falls  into  this  class.  The  other 
category  consists  of  small-scale  turbulence  (i.e.,  of  scale  length 
less  than  about  1000  feet).  Since  this  type  of  turbulence  involves 
motion  of  masses  of  air  which  are  roughly  the  size  of  the  sensor 
line  or  smaller,  its  effect  on  the  sensor  line  varies  from  one 
sensor  to  another.  Hence,  the  two  categories  are  distinguished  by 
whether  they  affect  all  the  sensors  simultaneously  in  the  same 
way.  The  final  component  of  sensor  output  is  the  inherent  measurement 
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error  in  the  instrumentation  itself.  Tvpically, it  consists  of  a bias 
effect  and  a random  component,  both  of  which  are  quite  small  if 
the  instrument  is  functioning  properly. 

The  effects  of  these  three  phenomena  are  additive,  and  they  result 
in  a pattern  of  sensor  outputs  which  is  typically  as  shown  in  Fig.  2.3. 

In  Fig.  2.3,  the  port  vortex,  located  at  approximately  +50  feet,  produces 
a minimum  in  the  sensor  located  at  +50  feet,  and  the  starboard  vortex, 
located  at  approximately  +325  feet,  produces  a maximum  in  the  sensor 
located  at  +300  feet.  Ambient  wind  biases  the  entire  pattern  by 
approximately  7 ft/sec,  and  small-scale  turbulence  produces  local  varia- 
tions in  the  pattern, such  as  the  local  maximum  at  -350  feet  and  the 
local  minimum  at  -200  feet. 

After  extensive  analysis  of  actual  sensor  tracking  data, it  has  been 
found  that  the  velocities  induced  by  vortices  are  only  measurable 
over  a region  of  about  150  to  200  feet.  In  other  words, the  effect  of  a 
vortex  appears  in  no  more  than  three  or  four  sensors  at  a time?  For 
example,  in  Fig.  2.3  the  starboard  vortex  only  affects  the  sensors  at 
250,  300,  and  350  feet.  Hence,  all  direct  information  about  vortex 
location,  at  any  particular  time,  is  embodied  in  only  a few  of  the 
instruments.  The  remaining  instruments  measure  the  combined  effect 
of  large-  and  small-scale  turbulence  and  ambient  wind. 

At  periodic  points  in  time,  all  sensor  outputs  are  sampled 
almost  simultaneously.  This  collection  of  20  or  so  sensor  outputs  is 
preprocessed  to  separate  the  various  etfects  of  large-  and  small-scale 
turbulence  from  the  effects  of  the  vortices.  The  vortex  effects 
manifest  themselves  as  a local  maximum  and  a local  minimum,  each  of 
which  is  concentrated  over  two  or  three  sensors.  The  regions  of 
vortex  influence  are  inferred  by  summing  the  outputs  of  pairs  of 
adjacent  sensors.  The  location  of  the  pair  of  sensors  which  indicates 
the  largest  sum  is  used  to  infer  the  region  of  the  starboard  vortex. 
Similarly,  the  location  of  the  pair  of  sensors  which  indicates  the 
smallest  sum  is  used  to  infer  the  region  of  the  port  vortex.  The 
outputs  of  the  sensors  in  each  of  the  two  pairs  are  ti.en  examined 
to  choose  a third  sensor  to  accompany  each  pair.  One  sensor  is 
chosen  so  that  three  sensor  outputs  associated  with  the  starboard  vor- 
tex form  a convex  group,  and  another  sensor  is  associated  with  the 
port  vortex  to  form  a concave  group.  For  example,  in  Fig.  2.3, the 

•Clearly,  this  situation  will  be  different  for  other  sensor  spacings. 
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FIGURE  2.3.  TYPICAL  SENSOR-INDICATED  VELOCITY  PATTERN  CAUSED  BY  VORTICES,  AMBIENT  WIND,  AND 
SMALL-SCALE  TURBULENCE  AS  A FUNCTION  OF  SENSOR  LOCATION 


sensors  at  ♦■300  and  +350  feet  form  the  maximum  sum  of  two,  so  that  the 
starboard  vortex  is  assumed  to  be  nearest  this  pair.  The  sensor 
at  +250  feet  is  appended  to  these  to  form  a convex  group  of  three. 

The  starboard  vortex  is  then  assumed  to  be  somewhere  between  +250 
and  +350  feet.  Simi lar ly. the  concave  group  of  sensors  at  0,  +50» and 
+100  feet  is  chosen, and  the  port  vortex  is  assumed  to  be  somewhere 
between  0 and  100  feet . 

To  determine  the  ambient  wind,  sensor  outputs  are  averaged. 
However,  the  average  does  not  include  data  from  the  group  of  three 
convex  and  the  group  of  three  concave  sensors.  In  other  words,  only 
th  se  5e."iSors  which  are  not  in  the  region  of  the  port  and  starboard 
vortice.'i  are  included  in  the  average.  Since  approximately  14  sensors 
remain,  a relatively  accurate  estimate  of  ambient  wind  is  thus  obtained. 
In  effect,  the  average  does  not  include  vortex  effects.  it  extracts 
the  large-scale  or  ambient  wind  effect  and  tends  to  suppress  the 
effects  of  small-scale  turbulence.  The  result  is  used  by  the  estimator 
as  a measure  of  the  ambient  wind  velocity  at  ground  level  (i.e.,  the 
variable  u in  Eg.  2.4). 

The  group  of  three  convex  and  the  group  of  three  concave  sensor 
outputs  are  then  used  to  infer  measurements  of  the  starboard-  and 
po’'t-vortex  positions,  respectively . First,  the  ambient-wind  velocity, 
inferred  from  the  averaging  process  described  above,  is  subtracted 
from  the  outputs  of  the  convex  and  concave  groups  of  three,  so  as  to 
eliminate  ambient  wind  from  the  data.  Then, the  groups  are  processed 
to  infer  measured  locations  of  the  port  and  starboard  vortices.  This 
processing  is  based  on  the  vortex  induced  velocity  model  of  Section 
2.3  and  explained  in  detail  in  Appendix  A.  The  formula  is 

X 


where 


(v  di  /v  -V  \ + V d^  /v  -V  \ + V d,/v  -v  \| 

I m^  1 \ >"2  '"3/  m2  2 ^ m^  m^^  j m^  3^  m^^  ”'2/ 

2 V d,/v  -V  \+  V d./v  -V  \+  V d,/v  -v  \] 

m ^ i y rn  2 ^ ^ f ^ 2 ~ \ ^ ^ ^ f ^ ^ *^  \ ^ ^ ^ 2 J 


, (2.7) 


x'  = inferred  measurement  of  vortex  position, 

'^m  * velocity  output  of  i*"^  sensor  of  the  group  of 

^ three  sensors  associated  with  the  vortex,  and 


dj^  = location  of  the  i 


. th 


sensor . 
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Equation  (2.7)  represents  the  lateral  position  of  the  vortex  an 
inferred  from  the  outputs  of  the  three  anemometers  wliich  are  lilcely 
bo  nearest  the  vortex.  As  described  above,  a number  of  factors 
corrupt  the  measurements  and  hence  produce  errors  in  inferred  position. 
As  a result,  the  inferred  measurements  are  processed  by  vor tex-tracking 
estimators  which  use  a history  of  measurements  to  provide  the  best 
estimates  of  vortex  position. 


2.5  VORTEX-TRACKING  ESTIMATOR 


As  discussed  in  Section  2.2  above,  the  vortex  motion  is  assumed 

to  bo  at  almost  constant  velocity  and  characterized  by  Eq.  2.4 

and  2.5.  The  horizontal  distance  x and  transport  velocity  v constitute 

the  two  state  variables  for  implementation  of  the  estimator.  Sensor  outputs 

are  sampled  periodically  and  processed  to  update  the  estimates  of 

these  two  state  variables.  Between  measurements,  the  vortex  motion  is 

extrapolated  in  a manner  which  is  consistent  with  the  state-variable 

equations  2.4  and  2.5.  Let  x , and  v , be  the  estimates  of  x and 

n-1  n-1 

V for  one  of  the  vortices  at  the  time  step  n-1,  based  on  measurements 

up  to  and  including  the  measurement  at  time  n-1.  Let  .x  and  be 

n n 

the  estimates  of  x and  v at  ,;he  time  step  n,  based  on  measurements  up 
to  and  including  the  measurement  at  n-1.  Then,  according  to  Eqs.  2.4  and  2.5 
and  the  fact  that  w is  a white  process  noise,  the  extrapolation  from 
time  n-1  to  time  n is 


X = 

n 


X , + 

n-1 


(2.8) 


(2.9) 


where  At  is  the  time  interval  between  n-1  and  n,  and  u^  is  the  ambient 
ground-wind  velocity  determined  from  the  ground-wind  sensors,  as 
described  in  Section  2.4  above. 


At  time  n,  a new  set  of  sensor  outputs  is  sampled.  As  described 
in  Section  2.4,  the  measured  positions  of  the  vortices  are  inferred 
from  these  sensor  outputs.  Let  j: ' be  the  inferred  measurement  of  the 
lateral  vortex  position  at  time  n.  Assuming  that  the  errors  in  successive 
inferred  measurements  are  independent,  the  Kalman -filter  theory  is 
applied  to  this  problem  to  yield  the  following  update  equations: 


X + K (x ' - X ) , 

n x„  n n ' 
n 


(2.10) 


I 
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V + 

n 


(x’ 

n 


(2.11) 


In  effect,  the  extrapolated  position  and  velocity  estimates  and 

V are  corrected  by  amounts  which  are  proportional  to  the  difference 
n 

between  the  measurement  and  the  extrapolated  position  estimate  x^ . 

The  proportionality  coefficients  or  filter  gains  and 

determine  the  character  of  the  filter.  According  to  the  Kalman 
filter  theory#  these  coefficients  are  determined  by  the  statistics  of 
the  measurement  errors  and  estimation  errors.  If  the  statistics 
indicate  that  the  measurement  is  much  more  accurate  than  the  current 
estimates,  then  the  coefficients  are  large,  giving  more  weight  to 
the  current  measurement  than  to  previous  data.  If,  however,  the 
statistics  indicate  that  the  measurement  is  much  less  accurate  than 
the  current  estimate,  then  the  coefficients  are  small  and  far  less 
weight  is  given  to  the  current  measurement  compared  to  previous  data. 

The  estimator  given  in  Eqs.  2.10  and  2.11  represents  a second- 
order  sampled-data  system  for  processing  measurement  data.  Application 
of  the  Kalman  filter  theory  provides  a direct  computational  method  for 
determining  the  filter  gains,  based  on  the  statistics  of  the  measurement 
errors,  the  strength  of  the  white  process  noise  w,  and  the  a priori 
statistics  of  the  initial  position  and  velocity  of  the  vortex.  For  the 
problem  at  hand,  none  of  these  statistics  are  well  known.  Furthermore, 
it  is  likely  that  they  change  with  time, so  that  to  apply  the  filter 
theory  directly,  one  must  provide  some  means  to  infer  these  statistics 
and  adjust  the  gains  accordingly.  In  addition,  it  is  important  to 
realize  that  the  basic  system  model  is  much  simplified, and, hence,  that  there 
are  inherent  modeling  errors. 

There  are  a number  of  possible  approaches  which  can  be  taken 
to  determine  appropriate  filter  gains.  The  ultimate  goal  is  to  choose 
gains  which  yield  good  system  performance  over  a wide  range  of  condi- 
tions. The  term  robustness  has  been  used  to  characterize  such  systems 
because  they  are  relatively  insensitive  to  changing  conditions.  To 
calculate  the  gains  by  applying  the  Kalman-f ilter  theory,  one  must 
select  covariance  parameters,  as  described  above,  and  solve  a matrix 
Riccati  eouation , the  solution  of  which  can  then  be  used  directly  to 
obtain  filter  gains.  The  gains  are  thus  determined  by  the  statistics 
assumed  when  developing  the  Ricatti  equation.  Experience  with  many 
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applications  of  filter  theory  to  real  problems  has  shown  that  typically, 
for  this  type  of  filterinq  problem,  the  qains  produce  a well  dampied 
filter  with  a frequency  response  which  balances  the  effects  of  measure- 
ment noise  aqainst  initial  uncertainty  and  the  effects  of  process 
noise.  With  this  fact  known,  it  has  been  determined  that  the  process  of 
quessinq  statistics  and  solvinq  the  Riccati  equation  is  a rather 
awkward  indirect  means  of  desiqninq  the  filter.  As  shown  in  Appendix  C, 
the  dampinq  ratio  of  the  filter  is  always  0.707.  Hence,  it  has  been  decided 
to  choose  and  so  that  this  dampinq  ratio  is  always  attained 

while  allowinq  experimentation  with  the  filter  bandpass.  This  parameter 
IS  then  determined  by  a cut-and-try  process,  usinq  actual  sensor  data 
to  yield  a practical  robust  system  which  attains  the  best  possible 
trackinq  performance.  Thus,  the  basic  filter  configuration  is  the  form 
obtained  from  the  Kalman  theory,  with  filter  qains  obtained  by 
experimentation  with  actual  vortex -sensor  data. 

Separate  filters  are  employed  for  each  of  the  two  vortices,  and  each  uses 
the  inferred  position  measurement  appropriate  for  the  vortex  it  is 
trackinq.  For  example,  the  starboard -vortex  tracking  filter  uses  the 
inferred  position  measurement  generated  from  the  group  of  three  convex 
sensor  measurements;  similarly,  the  port-vortex  tracking  filter 
uses  the  group  of  three  concave  sensor  measurements. 

The  vortex-tracking  filters,  as  described  above,  have  been  implemented 
in  a digital  computer  program  and  used  to  process  data  from  actual 
ground-wind  sensors  installed  at  London's  Heathrow  International  Airport. 
Appropriate  gains  have  been  chosen,  and  extensive  testing  has  been  done  to 
determine  the  effectiveness  of  the  trackers.  In  general  it  has  been  found 
that  in  relatively  calm  conditions,  when  winds  are  below  5 ft/sec,  it 
is  possible  to  choose  gains  which  make  the  trackers  work  quite  well. 

However,  when  there  is  appreciable  wind,  the  trackers  experience 
difficulty  in  initialization.  It  is  apparent  that  they  do  not 
adapt  to  the  wide  variations  in  the  locations  of  the  vortices  at  the 
times  when  useful  sensor  data  become  available. 

2.6  VORTEX-TRACKER  INITIALIZATION 

The  trackinq  problem  begins  when  an  aircraft  passes  over  the  sensor 
line.  Pressure  and/or  acoustic  sensors  are  triggered  by  the  aircraft  presence. 
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and  these  siqnals  alert  the  vortex -track inq  system.  After  the  air- 
craft passes  over  the  sensor  line,  the  vortices  roll  up  and  beqin 
to  descend.  There  is  often  considerable  variation  in  the  time  it 
takes  for  the  vortices  to  descend  to  the  point  where  useful  sensor 
data  become  available.  Aircraft  heiqht,  vortex  strenqth,  and 
meteoroloqical  conditions  are  all  factors  affectinq  descent  time. 
Furthermore,  the  wind  variations  with  altitude  provide  an 
additional  measure  of  uncertainty  as  to  vortex  location  when  they 
descend . 

A number  of  approaches  have  been  attempted  in  the  process  of  evolvinq 
a suitable  init ial ization  procedure,  including  adapting  the 
initial  estimates  to  the  ambient  wind,  cha.iging  filter  gains,  and 
delaying  the  tracker  initialization.  None  of  these  are  totally 
satisfactory.  Basically,  the  problem  is  one  of  detecting  the  onset  of  the 
vortex -induced  velocity  in  the  sensor  measurements,  and  initializing 
the  vortex-state  variables  (x  and  v)  to  values,  so  that  the  tracking 
estimator  can  lock  onto  the  vortices.  In  effect,  the  onset  of  the 
starboard-vortex-induced  maximum  and  port- vortex- induced  minimum, 
as  described  in  Section  2.4  and  illustrated  in  Figure  2.3,  must  be 
detected, and  the  tracking  estimators  initialized  accordingly.  Hence, 
it  has  been  determined  that  a measure  of  the  size  of  the  absolute  minimum 
and  absolute  minimum,  relative  to  the  average  sizes  of  local  maxima 
and  minima  in  the  other  sensors  in  the  line  would  be  useful.  To 
provide  this  comparison, a sample  standard  deviation  of  noise  in  the 
sensor  outputs  is  calculated  at  each  sample  time,  by  com- 
puting the  sum  of  squares  of  sensor  outputs,  subtracting  the  square 
of  the  sum,  dividing  by  the  number  of  sensor  outputs  used  and  taking 
the  square  root.  In  this  calculation  the  outputs  of  the  sensors 
comprising  the  maximum  sensor  pair  and  minimum  sensor  pair  as  described 
in  Section  2.4,  are  excluded.  The  resulting  sample  standard  deviation 
is  then  processed  through  a first  order  low  pass  filter  with  a time 
constant  of  6 seconds.  The  maximum  sensor  pair  and  the  minimum  sensor  pair 
are  then  each  divided  by  two,  the  average  wind  is  subtracted,  and  the 
results  are  passed  through  filters  which  are  identical  to  the  one 
employed  with  the  saunple  standard  deviation.  The  resulting  filtered 
averaged  maximum  sensor  pair  is  then  divided  by  the  filtered  sample 
standard  deviation  to  produce  a signal-to-noise  ratio  for  the  sensor 
data  associated  with  the  starboard  vortex.  Similarly,  the  filtered. 
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averaged  minimum  sum  of  two  is  divided  by  the  filtered  sample  stan- 
dard deviation  to  produce  a signal “to-noise  ratio  for  the  sensor  data 
associated  with  the  port  vortex.  These  two  signal -to-noise  ratios  are 
then  used  to  initialize  the  vortex-tracking  filters. 

Experiments  with  actual  sensor  data  determined  that  the 

signal~to-noise  ratios  must  be  greater  than  2.0  before  useful  tracking 

data  become  available.  Furthermore,  a delay  of  ten  seconds  after 

aircraft  detection  is  imposed,  during  which  time  no  initializations 

are  made,  to  allow  the  startup  transients  of  the  low-pass  filter 

to  settle.  Hence» no  track  initiations  are  done  until  these  conditions 

are  satisfied.  However,  once  the  delay  is  exceeded  and  a signal-to- 

noise  ratio  exceeds  two,  the  associated  filter  position-state  estimate 

X is  initialized  to  the  inferred  position  measurement  x',  and  the 
n n 

velocity-state  estimate  v is  set  to  zero.  In  effect, the  filter 
. ignores  all  past  data, and  places  total  faith  in  the  current  inferred 

position  measurement,  under  the  assumption  that  current  data  are  far 
more  important  than  previous  data.  The  signal- to-noise  ratios 
are  then  calculated  at  each  succeeding  measurement  time,  and  the 
I ■ changes  in  signal-to-noise  ratio,  from  one  sample  time  to  the  next, 

I are  monitored.  At  each  point  when  the  new  change  in  signal-to- 

! noise  ratio  exceeds  the  largest  change  observed  previously,  the 

1 estimator  is  re-initialized.  By  this  means,  the  rate  of  change 

signal-to-noise  ratio  is  monitored, and  the  estimator  is  initialized 
repeatedly  until  the  rate  of  change  of  the  signal-to-noise  ratio  begins 
' to  decrease.  This  procedure  is  terminated  at  40  seconds  after 

aircraft  detection,  under  the  assumption  that  the  vortices  should 
have  descended  and  produced  significant  sensor  signals  by  that  time. 

As  stated  previously,  the  track-initialization  procedure  evolved 
from  extensive  experimentation  with  actual  sensor  data.  It  has  been  found 
that  typically  the  signal-to-noise  ratio  builds  up  rapidly  as  the 
vortex  descends, and  then  decreases  gradually  as  the  vortex  dissipates. 
Monitoring  the  rate  of  change  of  signal-to-noise  ratio  is  found 
be  a reliable  moans  to  determine  the  onset  of  useful  tracking  data 
and  an  effective  indicator  of  when  to  re-initialize  the  estimators. 
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pf  2.7  BAD-DATA  SUPPRESSION  AMD  TRACKING-QUALITY  EVALUATION 

Experience  with  actual  tracking  data  also  has  eliminated  a problem 
with  sensor  measurement  errors.  At  times  when  the  siqnal-to-noise  ratio 
is  diminishing,  random  disturbances  in  the  sensor  outputs  can  produce 
very  large  variations  in  the  inferred  position  measurement  produced  by 
the  measurement  preprocessing  as  described  in  Section  2.4.  This 
problem  is  the  result  of  the  basically  nonlinear  relationship  between 
vortex  position  and  velocity  measurements.  It  has  been  found  that  these 
large  inferred  measurement  errors  produce  large  perturbations 
in  the  vortex  tracker,  resulting  in  possible  loss  of  the  true  vortex 
track  by  the  estimator. 

A straightforward  and  very  effective  remedy  has  been  found  by  simply 
limiting  the  size  of  the  corrections  which  can  be  made  to  the  esti- 
mate. If,  at  the  times  when  measurements  are  taken  and  preprocessed, 
the  inferred  measurement  of  vortex  position  differs  from  the  extra- 
polated estimate  of  position  by  more  than  200  feet,  then  the  measurement 
is  ignored  and  the  estimate  extrapolated  forward  without  a measurement 
incorporation.  By  this  device, a comparison  is  made  between  old  data 
and  new  data, and  if  there  is  a large  discrepency  the  new  data  are 
ignored.  Extensive  experimentation  has  shown  that  this  technique,  in 
concert  with  the  track -initiation  procedure  described  above,  is  most 
effective  in  suppressing  bad  data  and  seldom  ignores  good  data. 

A similar  procedure  is  used  when  the  vortices  approach  the 
boundaries  of  the  sensor  line.  When  this  occurs,  the  lack  of  sensor 
data  beyond  the  boundary  can  produce  large  errors  in  the  sensor  pre- 
processing procedure.  Hence,  once  the  estimated  vortex  track 
crosses  a sensor-line  boundary,  the  track  is  terminated. 

The  difference  between  the  inferred  measured  position  and  the 
extrapolated  estimate  of  position,  which  is  often  called  the  measure- 
ment residual,  has  also  been  found  to  be  an  effective  indicator  of  the 
accuracy  or  quality  of  vortex  tracking.  As  a means  of  quantifying 
this  indicator  the  residual  is  squared  at  each  measurement  incor- 
poration and  processed  through  a low-pass  filter.  This  procedure 
yields  an  approximate  running  evaluation  of  the  sample  mean  square 
value  of  the  residuals.  The  low-pass-filter  time  constant,  set  at  6 

seconds,  is  long  enough  to  provide  a measure  of  the  mean,  yet  short  i 


enouqh  to  adapt  readily  to  chanqing  tracking  conditions.  The  resulting 
tracking-quality  statistic  is  compared  against  fixed  threshold  values 
to  yield  an  evaluation  of  tracking  quality.  Letter  grades  are  given 
which  correspond  to  the  magnitude  of  the  square  root  of  the  statistic. 

The  thresholds  are  set  as  follows: 

A implies  O'-Zs”^  25  ft. 

B implies  25</s^  ^ 50  ft, 

C implies  50</s  < 75  ft, 

D implies  75 < /s  < 100  ft, 

E implies  100</s~^  150  ft,  and 
F implies  150</i~  , 

where  s is  the  filtered  mean-square  tracking-quality  statistic  calcu- 
lated as  described  above.  Thus,  a grade  of  A implies  excellent  tracking 
in  the  range  of  zero  to  25  feet  of  error,  and  a grade  of  F implies  very 
poor  tracking  with  deviations  larger  than  150  feet. 

2.8  TRACK  TERMINATION 

The  signal-to-noise  ratio  and  tracking-quality  indicators  are 
used  to  terminate  vortex  tracking.  As  the  vortices  age,  their  strength 
diminishes,  and  it  is  important  to  provide  some  definite  indication 
of  the  end  of  useful  tracking  data.  Experience  has  shovm  that 
typical ly,  after  the  40-second  delay  allowed  for  tracker  initializa- 
tions, the  signal-to-noise  ratio  generally  decreases.  When 
the  ratio  falls  below  two  and  more  than  40  seconds  have  elapsed  since 
aircraft  detection,  the  vortex  track  is  terminated.  Similar ly, when 
the  tracking-quality  indicator  falls  below  D,  so  that  the  rms 
residuals  exceed  100  feet,  and  more  than  40  seconds  have  elapsed 
since  aircraft  detection,  then  the  vortex  track  is  also  terminated. 

These  two  criteria,  by  which  track  terminations  are  made,  have  been  found 
to  be  effective  in  promptly  and  accurately  terminating  vortex  tracks 
when  useful  sensor  data  are  no  longer  available. 

Vortex  tracks  are  also  terminated  when  the  estimated  vortex 
position  passes  beyond  one  of  the  boundaries  of  the  sensor  line.  In 
this  case, the  rationale  is  that  useful  tracking  data  are  no  longer 
available  once  the  vortex  passes  outside  the  line  of  sensors. 


I 
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3.  VORTEX  TRACKING  RESULTS 


The  DOT/Transportation  Systems  Center  has  installed  ground-wind 
sensor  systems  at  a number  of  airports.  The  vortex  tracker  was  exten- 
sively tested  with  actual  data  tapes  supplied  by  TSC , 
containing  outputs  from  sensor  lines  deployed  at  London's  Heathrow 
International  Airport.  TSC  also  provided  vortex  track  printouts,  from 
its  own  computer  program,  to  be  used  as  a reference  and  a 
standard  of  performance  for  the  tracker  developed  at  CSDL.  Most  of 
the  testing  was  done  as  an  integral  part  of  the  process  of  tracker 
development.  Data  taken  on  three  separate  days  were  used  to  test 
the  system  over  a variety  of  tracking  situations.  The  first  data  tape 
(HM-94)  is  representative  of  a calm  day  with  winds  generally  below 
10  ft/sec.  The  second  data  tape  (HM-95)  was  recorded  on  a day  when 
large  headwinds  were  present.  Finally,  the  third  data  taoe  fHfl-62) 
represents  a relatively  high  crosswind  condition,  when  vortices  are 
blown  consistently  and  rapidly  across  the  sensor  line. 

In  the  following  sections,  selected  representative  vortex  tracks 
from  thase  three  tapes  will  be  presented.  A wide  range  of  aircraft 
from  Boeing  747' s to  BAC-lll's  is  included.  The  tracking  results  are 
discussed  by  way  of  comparisons  between  outputs  from  the  vortex- 
tracking estimator,  as  described  above,  and  outputs  from  the  tracking 
program  developed  at  TSC.  The  TSC  program  indicates  the  locations  of 
sensors  whose  outputs  yield  the  largest  number  of  maxima  and  the 
largest  number  of  minima,  from  samples  of  all  sensors  taken  at  1/7-second 
intervals,  during  successive  two-second  periods.  Since  these  maxima 
and  minima  are  effectively  the  raw  data  on  which  the  vortex -tracking 
estimator  operates,  the  comparisons  provide  a good  indication  of  how 
effective  the  tracking  estimator  is  in  processing  the  raw  data. 

Two  lines  of  sensors  were  deployed  at  Heathrow  Airport.  The 
so-called  outer  baseline  was  located  at  2400  feet  from  runway  thresh- 
old and  it  usually  senses  the  vortices  first.  The  inner  baseline  was  located 
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closer,  at  1475  feet  from  the  threshold.  Aircraft  pass  over  the  inner 
baseline  about  six  seconds  after  passing  over  the  outer  baseline. 

Aircraft  altitude  is  considerably  lower  over  the  inner  baseline  than 
is  the  case  over  the  outer  baseline.  Hence,  data  from  inner  baseline 
sensors  are  generally  of  higher  quality,  and  more  effective  tracking 
can  usually  be  done  using  these  data.  Vortex  tracks  derived  from  both 
sensor  baselines  will  be  presented  in  the  following  sections. 

i.l  VORTEX  TRACK I NO  IN  CALM  AIR 

* 

In  calm  air  vortex  tracking  is  relatively  easy.  Good 
tracking  data  are  usually  available  within  15  seconds  after  the  air- 
craft passes  over  a sensor  line.  The  vortex  motion  is  dominated  by 
ground  effect,  which  produces  an  almost  constant -velocity  motion  of 
each  vortex.  Sensor  signal -to-noise  ratios  are  generally  high, and 
accurate  tracking  is  possible  for  long  periods  of  time. 

Figures  3.1a  and  3.1b  present  a comparison  of  the  outputs  from 
the  TSC  program  and  the  tracking  estimator.  Sensor  data  from 
Heathrow  Tape  HM-94  (Case  11)  provide  the  input  to  the  two  programs. 

The  vehicle  which  generated  the  vortices  was  a British  Aircraft  Corp. 

Trident  airplane.  The  figures  present  data  processed  from  sensors 
in  the  outer  baseline. 

Figure  3.1a  is  the  output  from  the  TSC  program.  These  data 
can  be  viewed  as  representative  of  the  raw  data  on  which  the  vortex- 
tracking  estimator  operates.  Time,  designated  at  10-second  intervals 
after  aircraft  detection,  runs  vertically  down  the  left  hand  edge 
of  the  figure.  Distance  from  runway  centerline  is  presented  at 
intervals  of  100  feet  across  the  top  of  the  figure.  The  vertical 
column  of  apostrophes  in  Figure  3.1a  represents  the  runway  centerline, 
and  the  columns  of  dots  represent  the  edges  of  a corridor  of  ±150  feet 
in  width,  indicating  the  so-called  protected  zone  for  following 

aircraft^^.  The  letter  S designates  the  estimated  location  of  the  i 

starboard  vortex  and  the  letter  P represents  the  estimated  location  of 
the  port  vortex.  In  Pig.  3.1a  the  + sign  next  to  a P indicates  that  the 
port  vortex  may  lie  somewhere  between  the  P and  +.  Similarly  a side- 
by-side  S and  * indicate  a starboard-vortex  location  somewhere  be- 
tween the  two. 

As  seen  in  Fig.  3. la, good  data  from  the  port  vortex  appear  at  ■ 

about  15  seconds  after  aircraft  passage  and  persist  until  about  80  | 
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seconds,  at  which  point  the  data  become  erratic,  indicating  an 
unreasonable  discontinuous  jump  in  port-vortex  location  to  the  right 
of  the  starboard  vortex.  Good  tracking  data  on  the  starboard  vortex 
appear  at  about  20  seconds  and  persist  until  about  95  seconds, 
at  which  point  the  starboard  data  also  become  erratic.  The  port 
vortex  stays  within  the  +150  ft.  region  during  the  entire  80  second 
period;  the  starboard  vortex  passes  beyond  the  +150  ft  boundary 
at  some  time  between  40  and  60  seconds. 

The  columns  of  letters  running  vertically  down  the  left-  and  right-hand 
edges  of  Fig.  3.1a  are  measures  of  the  confidence  that  can  be  placed 
in  the  indicated  location  of  the  vortices.  For  example,  the  letter 
A in  the  right-hand  column  indicates  that  the  sensor  at  the  location  of 
the  G,  indicated  in  the  figure  at  the  same  point  in  time,  produced  the 
largest  output  at  all  14  samples  taken  in  the  proceeding  two-second 
interval.  Similarly  the  letter  B indicates  that  the  maximum  occurred 
in  that  sensor  90-99%  of  the  time  in  the  preceding  interval  and 
similarly  at  10%  increments  to  F,  which  represents  50-59%.  When  fewer 
than  50%  of  the  maxima  occur  in  one  sensor,  the  £ is  not  printed  and 
a blank  appears  in  the  right  hand  column.  Similarly,  the  left-hand 
column  applies  to  the  minima  which  are  associated  with  the  port 
vortices,  designated  by  P in  the  figure.  All  letter  grades  are  as 
indicated  above,  however,  when  the  minima  occur  less  than  50%  of  the 
time  in  one  sensor, during  a two-second  interval,  a dash  is  printed 
in  the  left  hand  column. 

Figure  3.1b  presents  the  port  and  starboard  tracks  as  indicated 
by  the  CSDL  vortex -tracking  estimator.  The  estimator  calculates  essentially 
continuous  tracks  to  high  resolution,  but  the  figure  quantizes 
the  estimates  to  within  50  ft.  intervals.  In  Figure  3.1b  the  column 
of  I letters  represents  runway  center  line,  lateral  distances  are 
designated  in  100-ft.  intervals  across  the  top  of  the  figure  and  the 
vertical  columns  of  dots  indicate  the  +150-ft.  protection  region. 

Letters  S and  P indicate  estimated  vortex  positions.  The  letter  Q 
indicates  a port-vortex  estimator  initialization,  as  described  in 
Section  2.6,  and  T represents  a starboard  vortex  estimator  initializa- 
tion. Thus,  the  port  vortex  is  initialized  at  16,  18,  20,  22  and  36 
seconds,  indicating  successive  increases  in  the  rate  of  change  of 
signal-to-noise  ratio  within  the  corresponding  two-second  intervals. 

The  port-vortex  track  stays  steady  at  about  -50  ft.  until  75  seconds, 
when  it  moves  out  to  about  -100  ft.  The  track  is  terminated  at  80 
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seconds  due  to  poor  tracking  quality.  Tracking  quality,  as  described 
in  Section  2.7,  is  printed  in  the  column  running  down  the  left-hand  side  of 
the  figure  for  the  port-vortex  estimator,  and  down  the  right-hand  side  for 
the  starboard-vortex  estimator.  As  can  be  seen, the  port-vortex  quality 
falls  abruptly  from  A to  E at  82  seconds,  precipitating  track  termina- 
tion in  response  to  the  erratic  behavior  of  the  sensor  data.  In  the 
period  from  about  20  seconds  until  80  seconds  the  estimates  are  very 
good,  as  indicated  by  the  A and  B grades  attained  in  that  period.  The 
starboard-vortex  track  is  initialized  at  26,  28,  32,  and  34  seconds. 

The  track  passes  beyond  the  +150  ft  boundary  at  54  seconds,  as 
evidenced  by  the  change  from  S to  * , in  Fig.  3.1b,  to  indicate 
starboard -vortex  position  outside  the  il50-ft  region.  If  the  port 
vortex  had  passed  outside  the  ±150-ft  region  the  P would  have  changed  to 
+ indicating  a port  vortex  positioned  outside  the  +150-ft  region. 

The  starboard -vortex  track  terminates  at  98  seconds  when  the 

sensor  data  become  erratic.  The  starboard-track  estimator  produces 
a smooth  track  with  high-quality  tracking  from  40  seconds  to  90  seconds. 
Track  termination  is  slightly  late,  as  evidenced  by  the  rather  abrupt 
and  unrealistic  motion  of  the  estimate  from  94  to  98  seconds. 

The  inner  baseline  data  taken  for  the  same  period.  Tape  HM-94 
(Case  11),  are  presented  in  Fia.  3.2.  Figure  3.2a 

represents  the  TSC  program  output, and  Fig.  3.2b  shows  the  tracking- 
estimator  output.  Performance  of  the  tracking  estimator  is  similar 
to  what  was  observed  in  Fig.  3.1.  in  this  case,  however,  when  the 
tracking  data  begin  to  deteriorate  after  120  seconds,  the  estimator 
tracker  is  able  to  continue  to  track  through  the  bad  period  and  the 
starboard  track  is  not  terminated  until  170  seconds.  Port-vortex  data 
begin  to  deteriorate  after  140  seconds,  and  the  tracking  estimator 
terminates  the  port  track  at  about  150  seconds.  It  should  be  noted, 
however,  that  even  though  a vortex  is  reported  after  120  seconds,  it 
is  likely  that  its  strength  is  greatly  diminished. 

Figures  3.3a  and  3.3b  present  data  for  the  outer  baseline  sensors 
with  a small  transport  aircraft,  in  this  case  a Douglas  DC-9.  These 
data  are  from  Tape  HM-94  (Case  12) . Good  sensor  data  are  available  until 
about  110  to  120  seconds, when  the  tracker  estimator  terminates  the 
tracks.  Tracking  quality  is  quite  high  over  the  entire  tracking  interval. 

Figures  3.4a  and  3.4b  represent  the  same  case  but  for  these  figures 
the  data  were  obtained  from  the  inner  baseline  sensors.  The  starboard- 
vortex  track  in  Fig.  3.4ta  illustrates  the  ability  of  the  tracking 
estimator  to  filter  out  sensor  errors  and  produce  a smooth  estimate  in 
the  face  of  rather  erratic  sensor  data. 
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Fiqures  3.5  and  3.6  present  trackinq  data  for  a Boeing  737 
transport.  These  data  are  from  Tape  HM-94  (Case  16).  The  port  track  in 
Fig-  3 . 5h  terminates  early  in  response  to  erratic  sensor  data  at 
about  65  seconds,  as  can  be  seen  in  Fig.  3.5a.  Comparison  of 
Figs.  3.6a  and  3.6b  demonstrates  the  ability  of  the  tracker  estimator 
to  lock  onto  a new  vortex,  as  its  effect  appears  in  the  sensors,  when 
the  remnant  of  an  old  vortex  still  remains  in  the  sensor  line,  as 
evidenced  by  the  port-vortex  indications  in  Fig.  3.6a,  at  -250  ft, 
up  until  about  16  seconds. 

3.2  VORTEX  TRACKING  IN  STRONG  HEADWINDS 

In  strong  headwinds  the  tracking  task  is  far  more  difficult  than 
in  calm  air.  There  is  sizable  turbulence  which  in  turn  produces  signals 
in  the  sensor  outputs  which  mask  the  vortex  tracks.  Furthermore,  the 
vortices  dissipate  quickly  and  the  process  of  track  termination  is  more 
critical.  On  the  other  hand,  the  signal -to-noise  ratios  are  generally  low,  and 
the  tracking  estimator  must  be  made  rather  sensitive  to  detect  marginal 
signals  so  spurious  signals  must  not  be  identified  as  vortex  tracks. 

Hence, a fine  balance  must  be  struck  between  sensitivity  and  the  ability 
to  reject  spurious  signals. 

Figures  3.7a  and  3.7b  present  outer  baseline  tracking  data  from 
Heathrow  Tape  HM-95  (Case  8)  with  vortices  produced  by  a Boeing  727 
aircraft.  Figure  3.7a,  which  is  the  output  from  the  TSC  program, 
serves  to  indicate  the  poor  quality  of  the  sensor  tracking  data,  due 
primarily  to  headwinds  of  about  25  ft/sec.  No  coherent  starboard 
vortex  track  is  evident,  and  only  the  trace  of  a port  vortex  track 
appears  on  Figure  3.7a.  Figure  3.7b  presents  the  output  from  the 
vortex-tracking  estimator.  No  starboard  track  is  reported, but  the 
estimator  is  able  to  track  the  port  vortex  from  about  34  seconds  to 
60  seconds.  Figures  3.8a  and  3.8b  present  the  TSC  program  outputs 
and  tracking-estimator  outputs,  respectively,  using  the  inner  baseline 
sensor  data  for  the  same  case.  Here  the  sensor  data  are  better,  and 
both  port  and  starboard  vortices  are  reported  by  the  tracking  estimator. 

The  vortices  are  acquired  promptly  at  about  15  seconds  and  are  terminated 
appropriately  at  about  34  seconds.  Tracking  quality  is  generally 
marginal,  as  evidenced  by  the  grades  received  while  tracking. 
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FIGURE  I. 7a  TSC  PROGRAM  OUTPUT  IN  STRONG 
HEADWINDS  USING  OUTER  BASE- 
LINE SENSORS  WITH  BOEING  727 
VORTICES 


FIGURE  3.7b  TRACKING  ESTIMATOR  OUTPUT  IM 
STRONG  HEADWINDS  USING  OUTER 
BASELINE  SENSORS  WITH  BOEING  727 
VORTICES 
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FIGURE  3.8a  TSC  PROGRAM  OUTPUT  IN  STRONG 

HEADWINDS  USING  INNER  BASELINE 
SENSORS  WITH  BOEING  727  VORTICES 
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FIGURE  3.9a  TSC  PROGRAM  OUTPUT  IN  STRONG 
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FIGURE  3.10a  TSC  PROGRAM  OUTPUT  IN  STRONG 

HEAOMINDS  USING  INNER  BASELINE 
SENSORS  WITH  BOEING  747  VORTICES 


FIGURE  3.10b  TRACKING  ESTIMATOR  OUTPUT  IN 
STRONG  HEADWINDS  USING  INNER 
BASELINE  SENSORS  WITH  BOEING  747 
VORTICES 


Figures  3.9  and  3.10  ate  again  a high "headwind  case  from  Tape 
HM-95  (Case  6);  with  a large  aircraft,  in  this  case  a Boeing  747. 

Figures  3.9a  and  3.9b  represent  the  outer  baseline  sensor  data,  pro- 
cessed through  the  TSC  program  and  the  tracking  estimator,  respec- 
tively. As  seen  in  Fig.  3.9a,  the  starboard  vortex  passes  outside 
the  sensor  line  in  about  40  seconds.  After  that  time  the  sensor 
at  +650  ft.  continues  to  have  the  highest  output,  resulting  in  the 
column  of  starboard  indications  at  +650  ft.,  from  40  seconds  until 
74  seconds.  The  estimator  tracker  rejects  this  spuriojs  signal  by 
terminating  the  starboard  track,  due  to  poor  tracking  quality,  at 
24  seconds.  The  termination  is  slightly  premature  in  that  the  track 
might  have  been  prolonged  until  about  34  seconds.  The  port-vortex 
track  continues  until  86  seconds,  at  which  point  the  signal-to-noise 
ratio  falls  below  2 and  the  track  is  terminated. 

Figures  3.10a  and  3.10b  reveal  the  tracking  results  for  the 
inner  baseline  sensor  data  for  the  same  case  (HM-95,  Case  6) . Signals 
are  strong  enough,  with  this  large  aircraft,  to  produce  excellent 
tracking  of  the  p>ort  vortex  and  reasonably  good  tracking  of  the  star- 
board vortex.  Erratic  sensor  outputs  cause  loss  of  tracking  quality 
and  termination  of  both  vortex  tracks. 

3.3  VORTEX  TRACKING  IN  STRONG  CROSSWINDS 

As  was  the  case  with  strong  headwinds,  strong  crosswinds  cause 
difficulty  in  vortex  tracking.  The  two  most  important  factors  affecting 
trackinq  nerformance  are  turbulence  and  the  fact  that  the  wind  may  displace 
the  vortices  by  hundreds  of  feet  from  the  runway  centerline  before  they 
begin  to  produce  useful  sensor  data. 

Figures  3.11a  and  3.11b  are  the  outputs  frcrni  the  TSC  program 
and  the  tracking  estimator,  produced  from  outer  baseline  sensor  data 
using  Tape  HM-62  (Case  11) . A light  transport  aircraft,  in  this  case  a 
Caravelle,  produced  the  vortices.  Sensor  outputs  are  generally 
erratic,  as  can  be  seen  in  Fig.  3 . 11a, and, hence,  no  coherent  vortex 
tracks  are  reported  by  the  tracking  estimator. 

Figures  3.12a  and  3.12b  show  the  results  of  processing  inner 
baseline  data.  In  this  case  the  starboard  data  are  quite  erratic,  but 
a port -vortex  track  is  discernable  in  Fig.  3.12a  from  12  to  40 
seconds.  This  signal  is  tracked  by  the  tracking  estimator  from  20 
seconds  until  32  seconds,  at  which  point  the  track  is  terminated  due 


to  poor  tracking  quality. 
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1 


Fiqures  3.13a  and  3.13b  illustrate  tracking  under  similar  condi- 
tions for  a Boeing  737.  The  data  are  the  outer  baseline  sensor 
outputs  from  Tape  HM-62  (Case  3) . Erratic  sensor  outputs  prevent 
effective  tracking  of  the  starboard  vortex;  however,  the  trace  of  a 
port-vortex  track  can  be  seen  in  Fig.  3.13a  and  the  tracking  estimator 
barely  catches  the  track  in  Fig-  3.13b. 

Figures  3.14a  and  3.14b  show  the  inner  baseline  results  for  the 
same  case.  The  port  vortex  does  not  produce  a sufficiently  high  signal- 
to-noise  ratio  for  tracking,  but  the  starboard  vortex  is  tracked  from 
16  seconds  to  34  seconds.  Tracking  quality  is  relatively  poor  due  to 
erratic  sensor  data. 

Figures  3.15a  and  3.15b  show  results  for  an  A-300  Airbus.  The 
data  were  obtained  from  the  outer  baseline  sensors  with  Tape  HM-62 
(Case  1).  The  tracking  estimator  produces  only  short  tracks.  The  delay 
in  initialization  of  the  starboard  track  is  the  result  of  a low  signal- 
to-noise  ratio  until  30  seconds. 

Finally,  Figs.  3.16a  and  3.16b  show  results  for  the  same  case 
using  data  from  the  inner  baseline  sensors.  Large  signal-to-noise 
ratios  produce  high-quality  tracking  over  most  of  the  tracking  period. 
Abrupt  loss  of  quality  is  due  to  passage  of  the  vortices  beyond  the 
sensor  line. 

3.4  VORTEX  TRACKING  RESULTS  SUMMARY 

The  results  presented  are  typical  of  many  vortex  tracks  produced 
With  real  data  from  the  Heathrow  tapes.  Typically  the  tracking 
estimator,  in  its  current  configuration,  is  conservative  in  the  sense 
that  vortex  tracks  are  initiated  somewhat  later  and  terminated  somewhat 
earlier  than  might  be  expected  from  performing  visual  pattern  recogni- 
tion on  the  raw  data.  This  type  of  performance  was  purposely  designed 
into  the  system,  as  a means  of  preventing  the  tracker  from  reporting 
false  vortices  in  response  to  turbulence,  gusts,  or  other  error  sources. 
The  tracker  is  quite  robust,  performing  well  over  a wide  range  of 
atmospheric  conditions  and  aircraft  types,  and  with  data  from  two  separate 
sensor  lines.  No  external  data  on  aircraft  type  are  required  by  the 
tracker  and  no  meteorological  data,  other  than  the  ground-wind  sensor 
data,  are  used.  In  calm  air  the  tracking  errors  are  generally  of  the 
order  of  25  feet  (rms) . In  turbulent  conditions  these  errors  can  range 
up  to  150  feet  (rms) , depending  upon  the  severity  of  the  turbulence 
and  the  strength  of  the  vortices  . 
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Experiments  were  also  run  with  sensors  removed  from  the  data, 
thereby  simulating  such  loss  of  information  as  would  occur  when  sensors  have 
failed  and  have  been  taken  off  line.  Removal  of  one  sensor  has  almost  no 
effect  on  the  tracking.  Removal  of  two  adjacent  sensors  introduces  some 
inaccuracy , but  the  estimator  is  able  to  maintain  the  track  across  the 
missing  sensors.  The  loss  of  three  adjacent  sensors  usually  causes 
difficulty  when  vortices  approach  the  region  normally  covered  by  the 
three  sensors.  Hence,  a substantial  loss  of  information  can  be  sustained 
before  the  estimator  becomes  ineffective. 

The  results  presented  above  pertain  solely  to  the  performance  of 
a tracking  algorithm  which  systematically  analyzes  sensor  data.  It 
would  be  a gross  error  to  attempt  to  evaluate  the  hazard  potential 
to  following  aircraft  solely  from  these  results.  Vortex  strength, 
altitude,  and  the  location  of  the  following  aircraft  with  respect  to 
the  vortices  cannot  be  obtained  from  this  algorithm;  these  parameters 
being  the  primary  elements  in  evaluating  hazard.  The  indication  of 
vortices  persisting  for  a long  period  of  time  is  purely  the  consequence 
of  the  algorithm's  ability  to  process  very  low-level  signals  on  calm 
days . 
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4.  SENSOR-FAILURE  IDENTIFICATION 


The  ground-wind  anemometers  and  their  associated  electronics 
are  subject  to  occasional  failures.  Typical  failure  modes  include: 
physical  damage  to  anemometers,  physical  misalignment  of  anemometers, 
bearing  wear,  loss  of  power,  60-cycle  interference,  scale  factor  errors, 
bias  shifts,  excessive  noise  in  output  signals,  and  many  other  possible 
difficulties . 

Although  there  are  many  failure  modes,  the  resulting  effects  on 
sensor  output  can  be  categorized  in  terms  of  a few  general  types  of 
failure  phenomena.  For  example,  physical  damage,  misalignment,  loss  of 
power,  and  scale-factor  error  often  manifest  themselves  primarily  as 
biases.  Similarly,  bearing  failure  and  60-cycle  interference  appear 
as  noise  in  the  output.  Hence, the  bias  shift  and  excessive  noise 
effectively  describe  a wide  range  of  failures, 

The  erroneous  signals  produced  by  failures  are  processed  as 
input  data  by  the  vortex  tracking  estimator.  Depenling  upon  the 
severity  of  a failure,  the  resulting  effect  on  tracking  performance 
can  range  from  a slight  reduction  in  tracking  accuracy  to  complete 
destruction  of  tracking  capability.  In  order  to  protect  the  vortex 
tracking  estimator  from  ingesting  bad  data,  a sensor  failure  detection 
and  identification  program  was  developed.  This  program  effectively 
compares  the  sensor  outputs  with  each  other,  utilizing  the  redundancy 
inherently  available  in  multiple  sensors  to  find  failures.  The 
program  looks  specifically  for  bias  failures  and  random  noise  failures 
in  the  sensor  output  data.  When  a failure  is  detected  and  the  bad 
sensor  identified,  its  output  is  henceforth  ignored  and  the  remaining 
sensors  are  used  to  track  vortices.  Although  there  is  a resulting 
loss  in  information,  the  tracking  estimator  has  consistently  demon- 
strated its  ability  to  track  vortices  effectively,  even  given 
failures  in  two  adjacent  sensors. 
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SENSOR  BIAS-FAILURE  IDENTIFICATION 


4 . 1 

The  sensor  bias-failure  identification  system  depends  upon 
comparisons  of  filtered  sensor  outputs  with  each  other.  No  failure 
identification  is  attempted  while  vortices  are  producing  strong  signals 
in  the  sensors.  This  is  accomplished,  in  the  failure -identification 
system,  by  not  processing  sensor  data  during  the  60  seconds 
following  aircraft  detection.  The  delay  effectively  gives  enough 
time  for  the  vortices  to  dissipate  to  the  point  where  they  are  no 
longer  a factor  in  failure  detection. 

At  sample  times  other  than  during  the  60  second  delay 
after  aircraft  detection,  the  output  of  each  sensor  is  processed 
through  a first-order  low-pass  filter.  There  is  one  filter  associated 
with  each  sensor, and  the  filters  are  all  identical,  with  time 
constants  of  200  seconds.  The  filters  suppress  high-frequency  noise 
but  pass  sensor  biases.  After  all  new  sensor  outputs  are  processed, 
the  filtered  outputs  are  averaged  to  obtain  the  sample  mean  across 
the  sensor  line.  Then, the  filtered  output  of  each  individual  sensor 
is  compared  with  the  sample  mean,  and  if  deviations  of  more  than 
five  ft/sec  are  found,  the  sensor  with  the  largest  deviation 
is  designated  as  failed.  The  filtered  output  from  that  sensor  is 
then  eliminated  from  the  sample  mean  and  the  mean  recalculated.  The 
comparisons  and  elimination  of  the  sensor  with  the  largest  deviation 
greater  than  5 ft/sec  are  repeated  until  all  deviations  are  less 
than  5 ft/sec. 

During  the  60-8econd  delay  allowed  for  vortex  dissipation, 
no  sensor  outputs  are  processed  and  no  failure  identification  is 
attempted.  To  allow  long-term  analysis  of  sensor  outputs,  the  filter 
state  variables  are  simply  held  constant  over  the  delay  period, 
thereby  effectively  carrying  the  past  history  of  sensor  performance 
across  the  delay  period. 

The  detection  and  identification  system  for  sensor-bias  failure 
was  tested  extensively  with  data  from  the  Heathrow  sensor  lines. 
Specific  periods  were  chosen  to  provide  data  containing  actual  sensor 
failures.  The  system  identified  all  known  sensor  failures  that  had 
been  previously  identified  by  visual  pattern  recognition  of  raw  data, 
and  no  false  identifications  were  made.  Ih  addition,  the  sys'^am 


identified  several  sensors  which  showed  marginal  f-er formance , which  had 
not  been  previously  identified  by  visual  pattern  recognition. 


Experience  has  shown  that  the  most  difficult  failures  to  detect 
and  identity  are  the  marginal  failures  which  exceed  the  5 ft/sec 
failure  threshold  by  only  a small  amount.  Typically  it  takes  two  to 
three  filter  time  constants  to  identify  these  failures.  Larger 
failures  take  less  time,  depending  upon  the  magnitude  of  the  failure. 
Hence  there  is  a compensatory  effect  in  that  the  larger  the  failure, 
the  more  quickly  it  can  be  identified  and  hence  the  shorter  the 
period  during  which  bad  data  are  ingested  by  the  estimator. 

4.2  SENSOR  NOISE-FAIL'JRE  IDENTIFICATION 

A procedure  similar  to  the  bias-identif ication  process  is  used 
to  identify  excessive  random  or  high-frequency  noise  in  sensor 
outputs.  The  60-second  delay  after  the  aircraft  identification  is  again 
used  to  eliminate  vortex  effects.  The  noise  identification  procedure 
calculates  sample  variances  of  the  sensor  outputs  and  compares  them 
with  each  other  to  identify  failures. 

At  sample  times  other  than  during  the  60-second  delay 
after  the  aircraft  detection,  the  output  of  each  sensor  is  squared  and 
the  result  processed  through  a first-order  low-pass  filter.  There  is 
one  filter  associated  with  each  sensor, and  the  filters  are  all 
identical,  with  time  constants  of  200  seconds.  The  output  of  each 
filter  is  a measure  of  the  sample  mean-square  output  of  the  associated 
sensor.  The  sample  mean  for  each  sensor,  obtained  from  the  sensor 
bias-failure  calculation,  is  squared  and  subtracted  from  the  sample 
mean  square  for  that  sensor,  to  obtain  a sample  variance  for  each 
sensor.  The  sample  variances  are  then  averaged  to  attain  a sample 
variance  for  the  entire  sensor  line.  Comparisons  are  then  made 
between  the  individual  sample  variances  and  the  average  sample 

2 

variance.  The  sensor  with  the  largest  deviation  above  25  (ft/sec)  is 
designated  as  failed  and  its  output  eliminated.  The  average  is 
recalculated  and  comparisons  repeated.  Iterations  continue  until 
all  deviations  are  less  than  25  (ft/sec)  . 
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This  noise -failure  identification  procedure  was  also  tested 
extensively  with  real  data  from  the  Heathrow  sensor  lines.  Behavior 
similar  to  that  of  the  bias -failure  identification  system  was 
observed.  All  failures  were  correctly  identified,  no  false  alarms 
were  issued, and  a number  of  marginal  failures  which  had  previously 
gone  undetected  were  found. 
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5.  COMPUTER  REQUIREMENTS 


The  algorithms  developed  for  vortex  tracking  and  failure 
identification  can  be  readily  applied  to  a system  deployed  in  the 
field.  The  algorithms  lend  themselves  nicely  to  modularization  and 
thus  can  be  readily  implemented  in  either  a central  processor  or 
a modular  computation  system  composed  of  many  small,  interconnected 
processors.  In  order  to  define  the  requirements  for  such  a system, 
the  algorithms  were  analyzed  to  determine  the  most  useful  modular 
configuration.  The  following  sections  describe  the  various  modules 
and  indicate  the  storage  requirements  estimated  for  each  module. 

5.1  AIRCRAFT  DETECTION  PROCESSING  AND  INITIALIZATION  FOR  NEW 

AIRCRAFT  ARRIVAL 

This  module  processes  data  from  the  pressure  sensors  which  are 
used  to  indicate  aircraft  arrival.  Certain  logical  checks  are  used 
to  reject  false  alarms.  Once  it  is  determined  that  a new  aircraft 
is  present,  the  program  provides  initialization  of  variables  in 
preparation  for  tracking  the  new  vortices  generated  by  the  aircraft. 
The  module  is  estimated  to  require  IOOj^q  words  of  storage. 

5.2  DATA-CONSISTENCY  CHECKS  AND  SCALING 

This  module  does  preliminary  processing  of  input  data  in  order 
to  prevent  bad -data  ingestion  by  the  estimator.  It  checks  for  large 
sensor-data  deviations  caused  by  bursts  of  electrical  noise, 
determines  loss  of  reference  voltages  and  performs  any  other  diagnos- 
tic tasks  that  do  not  require  long-term  filtering  as  utilized  in  the 
sensor  failure  identification  system  described  in  Section  4.0.  The 
module  also  scales  sensor  data  in  a manner  appropriate  for  subsequent 
utilization  in  the  vortex  estimator.  Storage  requirements  are  esti- 
mated at  120 words. 
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5.3  CLOCK-UPDATE,  DATA-OUTPUT,  AND  EXECUTIVE  FUNCTIONS 

This  module  performs  various  executive  functions,  including 
timekeeping,  output  data  processing, and  overall  system  management. 

Its  storage  requirement  is  estimated  at  250^^^  words. 

5.4  DETERMINATION  OF  LATERAL  POSITION  OF  SENSORS  PRODUCING  MAXIMUM 
OUTPUTS 

The  maximum  sum  of  two  sensor  outputs,  as  described  in  Section 
2.4,  and  the  lateral  p>osition  of  the  maximum  are  determined  by  this 
module.  A major  portion  of  the  module  consists  of  logic  for 
bridging  the  gap  between  sensors  when  a failed  sensor  is  present. 
Storage  requirement  is  estimated  at  200j^q  words. 

5.5  DETERMINATION  OF  LATERAL  POSITION  OF  SENSORS  PRODUCING  MINIMUM 
OUTPUTS 

This  module  is  essentially  the  same  as  5.4  above,  except  that  it 
searches  for  minimum  sensor  outputs. 

5.6  WIND  COMPUTATION  AND  VORTEX  PROPAGATION 

The  estimation  of  ambient  wind  velocity  is  done  in  this  module. 
The  estimate  is  then  utilized  to  extrapolate  vortex  estimated  position 
between  sample  times.  Storage  requirements  for  this  module  are 
estimated  at  80^^  words. 

5.7  DETERMINATION  OF  MEASURED  STARBOARD-VORTEX  POSITION  AND 
ESTIMATOR  UPDATE 

The  outputs  of  sensors  producing  the  maximum  output,  as  deter- 
mined in  5.4,  are  used  to  calculate  the  measured  position  of  the  star- 
board vortex.  This  measured  position  is  then  processed  by  the  esti- 
mator to  produce  updated  starboard  vortex  position  and  velocity  esti- 
mates. Storage  requirements  for  this  module  are  estimated  at 
words . 

5.8  DETERMINATION  OF  MEASURED  PORT-VORTEX  POSITION  AND  ESTIMATOR 
UPDATE 

This  module  is  the  same  as  5.7  except  that  it  performs  the 
same  tasks  for  the  port  vortex. 
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5.9  FAILURE  IDENTIFICATION 


This  module  implements  the  sensor's  failure-identification 
routine,  as  described  in  Section  4.  Storage  requirements  are 
estimated  at  I^Oj^q  words. 

5.10  COMPUTER  WORD  LENGTH 

The  variable  requiring  the  largest  dynamic  range  is  the  vortex 
location.  Its  greatest  magnitude  can  be  about  500  ft  and  the  resolution 
should  be  to  within  about  10  ft.  Hence  a dynamic  range  of  50  is 
required,  implying  about  7 binary  bits  of  information  plus  a sign  bit. 
Hence, an  8-bit  processor  is  sufficient  to  handle  the  problem  if 
multiplications  are  carried  out  to  16  bits  and  truncated  to  8 bits. 

5.11  COMPUTER  SPEED 

Any  modern  processor  with  basic  operation  times  of  the  order  of 
10  usee  is  sufficiently  fast  to  handle  the  problem. 
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6.  SIMULATOR  PLAN 


f 


A portion  of  the  effort  expended  in  this  research  included  the 
development  of  a plan  for  implementing  a simulation  of  the  vortex- 
tracking  problem.  This  simulator  will  be  used  to  aid  in  the  planning 
of  vortex-tracking  systems  for  airports.  The  simulator  is  to  be 
very  general,  allowing  great  flexibility  in  simulating  a wide  range 
of  tracking  situations.  Its  design  is  based  on  a modular  approach  to 
developing  the  various  required  elements  of  the  system.  The  modules 
are  chosen  in  a manner  consistent  with  the  various  physical  phenomena 
or  operational  elements  that  influence  vortex  tracking.  Their  size 
is  based  on  the  desire  to  have  each  module  as  self-contained  as 
possible  with  a minimum  of  input  and  output  requirements.  In  the 
sections  that  follow,  the  roles  of  the  various  modules  will  be 
discussed  and  the  overall  simulator  system  design  explained.  A diagram 
of  the  simulator  configuration  is  presented  in  Figure  6-1. 

6 . 1 Aircraft  Module 

The  initial  locations  of  vortices  are  determined  by  the  path 
flown  by  the  aircraft  generating  the  vortices.  This  simulation  module 
generates  the  flight  paths  of  aircraft  on  approach  to  landing  and 
on  climb-out  after  takeoff.  It  effectively  provides  the  initial 
conditions  for  the  vortex  simulation  module,  which  is  described  in 
the  next  section. 

The  module  will  also  contain  a simplified  aircraft  dynamic 
simulation  along  with  a closed-loop  simulation  of  the  pilot  or  auto- 
matic landing  system.  The  purpose  here  is  to  allow  some  evaluation 
of  the  effect  of  a vortex  encounter  on  an  aircraft  during  final 
approach  or  climb-out. 
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6.2  VORTEX-DYNAMICS  MODULE 

This  module  will  produce  a simulation  of  vortex  motion  and 
determine  the  velocities  produced  by  the  vortices  at  designated 
points  in  space.  All  imp)ortant  aspects  of  vortex  dynamics  will  be 
included  so  as  to  provide  an  accurate  representation  of  the  motion. 

The  simulation  will  be  done  in  three  dimensional  space,  with  the 
locus  of  vortex  centers  represented  as  connected  line  segments. 

Vortex  motions  will  be  represented  by  motions  of  these  lines  in  the 
three-dimensional  space.  The  initial  location  and  shape  of  each 
vortex  line  is  determined  by  the  aircraft  simulation  program 
described  above.  Ground  effect,  ambient  wind,  air  turbulence,  vortex 
interactions,  and  vortex  decay  will  all  be  included  by  appropriate 
mathematical  modelling. 

Also  stored  within  this  module  will  be  the  geometry  of  the 
particular  airport  being  studied.  Provision  will  be  made  to  specify 
the  airport  terrain,  locations  of  runways, and  placement  of  vortex 
sensors.  With  this  data  and  the  lines  of  vortex  centers,  as 
described  above,  vortex  induced  velocity  at  desired  points  will  be 
generated.  In  particular, the  induced  velocities  at  the  locations  of  the 
vortex  sensors  and  at  nearby  runways  can  be  supplied  as  output  from 
the  program. 

6.3  SENSCR-SIMULATION  MODULE 

This  module  provides  mathematical  models  of  the  ground-wind 
anemometers.  Its  input  is  the  local  air  velocity  at  locations  of  each 
of  the  sensors, and  its  output  is  indicated  velocity,  including  sensor 
errors.  These  errors  will  include  a bias  term,  scale -factor  devia- 
tions and  a high~frequency  noise  component.  Provision  will  also  be 
made  to  simulate  misalignment  of  the  sensor  relative  to  its  desired 
orientation . 

Simulations  of  sensor  failures  will  also  be  implemented  in  this 
module.  A wide  range  of  possible  malfunctions  will  be  provided, 
in  order  to  allow  evaluation  of  tracking  effectiveness  in  the 
presence  of  various  failures.  The  type  of  failure,  magnitude,  and  time 
of  repair  can  all  be  specified  by  appropriate  Inputs  to  the  module. 
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TRACKINn-SYSTEM  MODULE 


In  an  actual  airport  installation  the  tracking  task  will  be 
accomplished  in  a dedicated  digital  computer  system.  It  is  currently 
envisioned  that  this  computer  system  will  be  of  modular  design,  made 
up  of  a number  of  microprocessors  tied  together  by  appropriate  data 
communication  buses. 

The  simulator  will  be  designed  to  allow  simulation  of  this  type 
of  system.  Both  hardware  and  software  will  be  implemented  as  elements 
within  the  simulator;  this  will  allow  evaluation  of  the  data-pro- 
cessing  system's  hardware  configuration  as  well  as  checkout  of  tracking 
algorithm  software.  By  providing  a detailed  simulation  of  the  data 
processing  system,  it  is  anticipated  that  sizable  savings  in  hardware 
and  software  costs  can  be  realized  by  appropriate  system  checkout  on 
the  simulator. 

6.5  DATA  OUTPUT-AND-DISPLAY  MODULE 

A very  carefully  designed  data  output-and-display  system  will 
be  utilized  to  provide  simulation  results  to  users.  The  system  will 
simulate  the  display  system  which  will  be  available  to  air-traffic 
controllers;  its  purpose  is  to  inform  the  controllers  of  whether 
there  are  any  vortices  within  the  + 150-ft  protected  region. 

This  simulator  module  will  also  provide  users  with  a versatile 
means  of  monitoring  variables  generated  within  the  simulator,  for 
purposes  of  evaluating  system  performance.  Users  will  be  able  to 
choose  from  a series  of  fixed  data  formats,  which  will  allow  a wide 
range  of  choice  relative  to  the  amount  of  detail  which  will  be  pro- 
vided in  the  output.  In  addition,  provision  will  be  made  to  allow 
users  to  construct  a desired  output  format,  which  might  differ  from 
one  of  the  fixed  formats.  A very  versatile  and  carefully  implemented 
system  for  control  of  data  output  will  be  devised. 

6.6  EXECUTIVE  PROGRAM 

The  various  simulation  modules  will  be  tied  together  through 
an  executive  program  which  will  serve  to  organize  the  modules,  exer- 
cize logical  decision  making,  and  provide  the  overall  management  of  the 
simulator  while  it  is  running. 


6-3 


New  Aircraft 


FIGURE  6.1  SIMULATOR  CONFIGURATION 


7.  CONCLUSIONS 


An  effective,  automatic  vortex-tracking  software  system  has  been 
developed  to  process  data  from  ground-wind  propeller  anemometers. 

System  performance  has  been  verified  using  actual  sensor  data  gathered 
at  London's  Heathrow  International  Airport.  The  system  works  well  over  large 
variations  in  wind  conditions  and  variations  in  strength  of  vortices. 

It  is  tolerant  of  sensor  errors  and  adapts  readily  to  loss  of  sensors 
due  to  malfunction. 

An  effective  software  system  for  detecting  and  identifying  sensor 
failures  has  also  been  developed.  Actual  sensor  data  were  used  to 
verify  performance  of  the  failure  detection  and  identification  software. 

All  failures  previously  identified  by  visual  pattern  recognition  with 
raw  data  were  also  recognized  as  failures  and  identified  by  the 
automatic  system,  and  no  false  alarms  were  given.  In  addition, certain 
marginal  failures,  which  had  previously  gone  undetected  by  visual 
pattern  recognition,  were  found  by  the  automatic  system. 

The  software  has  proven  itself  with  actual  sensor  data  under  a 
wide  range  of  operating  conditions.  It  can  be  readily  adapted  to 
real-time  tracking  and  for  application  to  future  vortex  warning  systems 
installed  at  airports.  Initial  steps  have  been  taken  to  modularize  the 
software,  in  preparation  for  its  implementation  in  either  a central 
data  processing  installation  or  a modular  computer  system. 

A tentative  plan  for  development  of  a general-purpose  vortex - 
tracking  simulation  program  has  also  been  presented.  This  simulator 
will  be  an  effective  tool  for  the  study  and  development  of  vortex 
warning  systems.  If  properly  designed,  the  simulator  will  be  useful 
for  planning  sensor  Installations,  software  development  and  verification, 
performance  evaluation,  and  for  studying  various  combinations  of  sensor 
deployment,  operational  procedures,  and  sensor- failure  detection. 
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APPENDIX  A 

DERIVATION  OF  EQUATION  FOR  INFERRED 
VORTEX  POSITION  MEASUREMENT 


A useful  equation  can  be  derived  for  the  lateral  position  of  a 
vortex  in  terms  of  three  measurements  of  induced  velocity.  The  equation 
is  based  on  an  assumed  i velocity  field  for  the  vortex  induced  velocity. 
From  Eq . 2.6  the  horizontal  velocity  measured  by  sensor  i is 


rh 


nlh^  + (x-d^)^l 


(A-1) 


I 

i 

i 


r = vortex  circulation 
h = vortex  altitude 
X = vortex  lateral  position 
d^  = sensor  lateral  position. 

If  there  are  three  sensors,  then  three  equations  of  the  form 
(A-1)  can  be  written  for  the  three  sensors  with  i = 1,  2,  3. 

r h 

The  equations  for  sensors  1 and  2 can  be  combined  to  eliminate 
obtaining 


-2d,x+d,^j=v  (h^+x^) -2d,x  + d 

1 1 m.,  2 


(A-2) 


and  similarly  for  sensors  1 and  3 
(h^+x^)-2dj^x  + d^ 


s:  V. 


(h^+x^)  - 


2d3X  + dj 


(A- 3) 


These  two  equations  can  then  be  combined  to  eliminate  the  term 
2 2 

(h  tx  ),  obtaining 


V -V 

">1  >"3 


2 2 

X /2d,v  -2d_v  \+  V d-,  - V d, 

I 1 mj^  2 I mj  2 mj^  1 

x/2d,v  -2d-V  W V d^  - V ? 
y 1 mj^  3 m^  3 mj^  1 


(A-4) 


A-1 


which  can  then  be  solved  for  x to  yield 


APPENDIX  B 


SUBROUTINE  SUMMARIES, 

PROGRAM  FLOWCHART,  AND 
PROGRAM  LISTINGS 

The  subroutine  summaries  are  preceded  by  a glossary.  The  basic 
format  of  each  subroutine  summary  appearing  on  the  following  pages 
is  composed  of  the  following  terms: 

A)  Purpose 

B)  Calling  Sequence 

C)  Inputs 

D)  Outputs 

E)  Subroutines  Called 

F)  Subroutine  Length 

G)  Warning  or  Error  Messages . 

A flowchart  is  then  enclosed  for  the  significant  and  complicated 
subroutines,  followed  by  the  FORTRAN  listing  for  each  subroutine. 


ESTIMATOR -TRACKER  GUDSSARY 


ABSPL  Absolute  value  of  port-vortex  position 

ABSPR  Absolute  value  of  starboard -vortex  position 

alt  Altitude  of  port  or  starboard  vortex 

Cy  Velocity-filter  measurement  weight 

Cjj  Position-filter  measurement  weight 

DELT  Time  step 

DENOM  Denominator  used  in  calculating  the  measured  lateral 

position  of  either  the  port  or  starboard  vortex 

DMN  Locations  of  the  three  smallest  average-velocity  sensors 

at  time  N 

DMX  Locations  of  the  three  largest  average-velocity  sensors 

at  time  N 

DN  Lateral  positions  of  the  three  sensors  which  represent 

the  minimum  average  velocity 

DSQ  Used  to  determine  if  the  magnitude  of  the  difference 

between  the  measured  position  and  the  estimated  position 
is  consistent  with  expected  measured  and  estimated  errors 

DX  Lateral  positions  of  the  three  sensors  which  represent 

the  maximum  average  velocity 

Dl  Lateral  position  of  the  sensor  with  the  largest  velocity 

for  starboard,  smallest  velocity  for  port 

D2  Lateral  position  of  the  sensor  with  the  smallest  velocity 

for  starboard,  largest  velocity  for  port 

I Index  pointer 

IBLMN  Indicates  that  the  minimum  value  for  the  sum  of  two  sensors 

is  on  the  left-hand  boundary 

IBLMX  Indicates  that  the  maximum  value  for  the  sum  of  two  sensors 

is  on  the  left-hand  boundary 
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0 


IBRMN 


IBRMX 

ICL{0,I}* 

ICLK 

ICB{0,I) 

IDX 

IFST 

IFSTI 

II 

IPDX 

IPX 

IS 

JJ 

J1 

J2 

J3 

J4 

L 

LARG 

MNL{0,I) 

MXR(0,I) 


Indicates  that  the  minimum  value  for  the  sum  of  two 
sensors  is  on  the  riqht-hand  boundary 

Indicates  that  the  maximum  value  for  the  sum  of  two  sensors 
is  on  the  right-hand  boundary 

Indicator  used  in  termination  logic  on  port  side 
Internal  clock  used  to  count  frames  per  second 
Indicator  used  in  termination  logic  on  starboard  side 
Index  pointer 

Indicator  used  on  outer  baseline  to  signal  arrival 
of  new  aircraft 

Indicator  used  on  inner  baseline  to  signal  arrival  of 
new  aircraft 

Counter 

Index  pointer 

Index  pointer 

Sensor  number 

Counter 

Sensor  location  on  the  left-hand  side  of  the  minimum  or 
maximum  sensor-pair  chosen 

Sensor  location  on  the  left-hand  side  of  the  minimum  or 
maximum  sensor-pair  chosen 

Sensor  location  on  the  right-hand  side  of  the  minimum  or 
maximum  sensor-pair  chosen 

Sensor  location  on  the  right-hand  side  of  the  minimum  or 
maximum  sensor-pair  chosen 

Measured  lateral  position  of  the  port  vortex 
Large  floating-point  number  used  for  Initialization 
The  location  of  the  first  sensor  on  a baseline 
The  location  of  the  last  sensor  on  a baseline 


0 - means  outer  baseline 

1 - means  inner  baseline 
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NP  Number  of  sensors  minus  the  number  of  failed  sensors 

minus  three 

R Measured  lateral  position  of  the  starboard  vortex 

RMSVWfO.l)  Root-mean  square  wind  velocity 

SCONST  Scaling  constant 

SECC  Forty-second  threshold 

SECO  Ten-second  threshold 

SEC60  Sixty-second  threshold 

SMALL  Small  floating-point  number  used  for  initialization 

SMAX  Maximum  velocity  for  the  sum  of  two  sensors 

SMAX3  Strength  of  the  starboard  vortex 

SMIN  Minimum  velocity  for  the  sum  of  two  sensors 

SMIN3  Strength  of  the  port  vortex 

SNLBST  Intermediate  variable  used  to  store  the  largest  difference  between 

two  successive  signal-to-noise  ratio  values  on  the  port  side 

SNL{o,l}  Signal  to  noise  ratio  for  the  port  vortex 

SNLP  Intermediate  variable  used  to  store  a previous  signal-to-noise 

ratio  value  on  the  port  side 

SNRBST  Intermediate  variable  i^sed  to  store  the  largest  difference 

between  two  successive  signal-to-noise  ratio  values  on  the 
starboard  side 

SNR{0,l}  Signal-to-noise  ratio  for  the  starboard  vortex 

SNRP  Intermediate  variable  used  to  store  a previous  signal-to-noise 

ratio  value  on  the  starboard  side 

SPRTL  Print  trigger  to  indicate  start  of  the  port-vortex  track 

SPRTR  Print  trigger  to  indicate  start  of  the  starboard -vortex  track 

SRMSVW  Smoothed  root -mean -square  wind  velocity 

SSMN  Smoothed  signal-to-noise  ratio  for  the  port  vortex 

SSMX  Smoothed  signal-to-noise  ratio  for  the  starboard  vortex 


0 - means  outer  baseline 

1 - means  inner  baseline 


SSQVW  Sum  of  the  squares  of  good  sensor  velocities,  excluding 

the  two  sets  of  sensors  with  the  largest  and  smallest 
velocity 

SUM  Intermediate  variable  used  to  store  the  maximum  or  minimum 

velocity  for  the  sum  of  two  sensors 

TEM  Intermediate  variable  used  in  the  calculation  of  the 

measured  lateral  position 

TEMP  Temporary  variable 

TPIE  Constant 

V Velocity  estimate  of  the  port  vortex 

Li 

V„  Velocity  estimate  of  the  starboard  vortex 

V^  wind  velocity  estimate 

WOLD  g-At/Time  Constant  smoothing  old  data 

WNEW  (1-WOLD)  used  for  smoothing  new  data 

X Scaled  vortex  motion  as  recorded  by  sensors  minus 

wind  estimate 

XHAT  Scaled  vortex  motion 
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ESTIMATOR-TRACKER  SUBROUTINE 


A.  Purposes 

To  track  automatically , and  to  record  measured  vortex  positions  as 
a function  of  time. 

B . Calling  Sequence 

CALL  EST  (Z,  ISPIKE,  NS,  ISN,  ISDST,  lORD,  NETS,  NF,  P, , P_, 

Li  R 

SSQL,  SSQR,  SMHL,  SMXR,  SPRTL,  SPRTR) 


Inputs 

Through  variables  in  calling  sequence 

Z^  Ground-wind  sensor  data  as  recorded  by  propeller 

anemometers 

ISPIKE  Trigger  to  indicate  spike  in  GWS  data 

NS  Total  number  of  sensors  in  a baseline 

Indicates  failed  sensors  in  a baseline 

Sensor  lateral  positions  in  a baseline  with  respect 
to  the  centerline  of  a runway 

Sensor  locations  in  a baseline  with  respect  to  the 
centerline  of  a runway,  ordered  from  port  to  starboard 

Ambient  meteorological  wind  data  recorded  at 
various  levels  on  multiple  towers 

Total  number  of  failed  sensors  in  a baseline 
The  remaining  inputs  are  passed  through  labeled  common  blocks. 
/INOU/  -Specifies  input,  output,  and  scatch  device  numbers. 

/INPUTl/,  /PASS21/,  /ORDERS/  -See  glossary  for  definition  of 
variables  in  these  common  blocks. 


ISN^ 

ISDST^ 


lORD, 


METS, 


NF 


) 
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D.  OUTPUTS 

Through  variables  in  calling  sequence 
P Position  estimate  of  the  port  vortex 

Li 

P^  Position  estimate  of  the  starboard  vortex 

SSQL  Smoothed  quality  indicator  for  the  port  vortex 

SSQR  Smoothed  quality  indicator  for  the  starboard  vortex 

SMNL  Smoothed  strength  of  the  port  vortex 

SMXR  Smoothed  strength  of  the  starboard  vortex 

SPRTL  Print  trigger  to  indicate  start  of  the  port  vortex  track 

SPRTR  Print  trigger  to  indicate  start  of  the  starboard  vortex  track 

The  remaining  outputs  are  passed  through  a labeled  common 

block. 

/EST{0,I}/  -See  glossary  for  definition  of  variables  in  this 
common  block. 


E . Subroutines  Called 

CK,  FDI 
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VORTEX  ESTIMATOR  TRACKER  MODULE 


I 

i 


1 
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Check  If  measured 
position  r minus  esti- 
mated position  Pr  <s 
greater  in  magnitude 
than  200  ft. 


Update  position  and 
I velocity  estimates 

P,  ■ f(P|,(n).C,.P) 
V,  ■ f(v„ln).C,.«) 


Logic  to  determine 
when  to  start  star- 
board-vortex track 


Logic  to  determine 
when  to  terminate 
starboard-vortex 
track 


Compute  the  measured  lateral 
position  (L)  of  the  port  vortex, 
the  quality,  the  altitude,  and 
the  strength,  using  the  group 
of  three  concave  sensors 


Check  whether  measured 
position  L minus  esti- 
mated position  P,  Is 
greater  in  magnitude 
than  200  ft. 
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^0 

I Update  position  and 
' velocity  estimates 



Logic  to  determine 
when  to  start 
port-vortex  track 

_ : ':i 

I Logic  to  determine 
I when  to  terminate 
, port-vortex  track 


Has  60  seconds  expired 
since  the  last  aircraft 
arrival 


CALL 

Failure-Detection  Module 


RETURN 


BtSl  AVAILABLE  COPY 


ESTIMATOR -TRACKER  SUBROUTINE 


SUbBT  IT  INE  FSTI  Z,  lEPlK^.NS. ISN. I SrST , ICRD .ME TSt NF tPL .PC t «SQL • SSUR 
I .SMNL .S  MHP . SPPtL. SPPtPt 
C=’*L  L.  LAP'..  MS,  Ml 

InTEG^p  DMN.  P'IX.  Z.  SitCC.  secc.  SPFTF.  SPRTU.  SEC60 

riMFNSJ  ON 

1 ONITI,  r>Xl-J).  OMNIA).  CMXIA).  tem|?),  XI!).  XHATOO) 

or  MEN?)  ON 

) lENl’O).  I'-03T|30).  I'^Pri?").  METSIIC).  ZI30) 

COMMON/  INQU/ICP  ,IPP  . IP'J.  lOT  . I E 1 . I E *> . I T V . I S3  . t S«  . 1 . I S6 

COMMON/ INTT) /Cx .CV.OSO.DELT.SSCO.EFCC  .MfLC.MNEW 
common /PASS 21 / I FST. IFSTl 
CCMMrN/CCOE02/MNLO.MXFO 

comm Ofl/ESTn/PMJVMO.  SNLD  .SNPC  .ICLK.  ICPC.  ICLO 
CATA  SEC<AC/4?0/ 

CATA  LAPG/1.E2S/ 

TATA  SMALL/-!.  - 2^/ 

CATA  TP |E/s,2A )l AJ. oe/ 

C COUIVAlFNCE  |MrTS|P).F3) 

SCnNS’»I0C.0/IAETSI9)-METSI’ 0) ) 

IF  I I •'ST  ,NE.  I ISO  TO  MO 

C initialization  for  every  NFM  aircraft  arrival 

CO  100  I«1.NS 
XHAT I I ) *C .0 
100  CONTINUE 
VV.O."' 

VRiO.'' 

VC»0.C 
pRxSO.O 
PL*-SC. P 
IFST*0 
ICLK*0 
SSQRcC. 0 
«SQL*0. 0 
Smnl*0. 0 
SMXB  »0.C 
SPMS  VW*  0. 0 
«SMX»0. 0 
SSMN.O. 0 
SNRaST*  0. 0 
SNLBSTa 0, 0 
3PPTPX  0 
S^RTl *0 
IFP0«A 
irL'jAA 

no  IF|  ISPIKE.NF.OIGO  TO  26  S 
C UPfOATE  PREP  iL  TCP 
rn  iz"  i«i.NS 
ICX*  inpci  I » 

IF| I SNI lOXJ.FO.l )Gn  TO  120 
XHAT II ) ■! FLOAT I Z|  tDX )-METS I A ) I ) ASCONST 
120  CONTINUE 

C FINC  max  and  min  MrAN  VELOCITY  MEASUREMENTS 
SMAXaSMALL 
SMtN>LARG 
iaPMX»0 
inc  mni-i? 
np«NS-3-nf 
1^1 
Ji  ■! 

125  lOx^IOPOIJI) 

1F|  ISNI  IDXl.EQ.DGO  TC  130 
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GO  TO  40 

130  I3«Jl4t 

CALL  CK ( I S. ISN. lanO.NSt Jl » 

UO  IOX<>IC<RD(  J14n 

I'=^(  I AM  irxt.cQ.  1»G0  TO  ISO 

GO  TO  1 60 
ISO  IS>JU2 

CALL  CK  ( I S«  t SN«  IORO«NA<  ) 

|60  I0X=lnRDlJ2*l  » 

IR  ( ISN(  IDXI  (EG.  1 )G0  TC  ! ''0 
J3»J2f 1 
CO  TO  1 60 
iro  IS=J2t2 

call  CK I I 3> I SM. IORD (NSt J3 > 

180  !OX«>!nRC(  J3TI  » 

1P< ISNC  1PX»,E0«1 »GO  TO  185 
J4*J3+1 
GO  TO  2 0C 
135  IS-J3+2 

call  CK ( I St ISNi lOROtNSf J4I 
200  SUM=XHAT< J2»TXHAT(J1J 

IF|  SUM.  LE  .S><AX)  GO  TC  205 
S><AX  = 5UM 
tRLMXcO 
OMXC 1 )s J1 
CMX<  2I  = J2 
DMX( 3 J3 
OMX( 4)3 J4 
IF( 1*EQ«1 IIBLHX^I 
IF( ! •FQ.NPI I BRMX«1 
205  !F(  SUM,  GF  ,S'^IN»  GO  TO  210 
SMIN=FUM 
tBLMN=0 
OMN( 1 )= Jl 
CMN(2)= J2 
OMN{ 3»» J3 
DMN(4 )= J« 

IF<  t ,60,1  IIBLMN^I 
IFI  I ,6Q,NPi IBRMNal 
210  IFI  1 .GE  ,NPlGn  TO  215 
♦! 

J!>jl4l 
GO  TO  1 2S 
215  JI-DMX(l) 

J2«DMX(  2) 

J3«DMX( 31 
JAaOMXI  4) 

IF|  XHAT( J3) ,LF, XHAT ( J2I IGO  TC  22C 
CMX( 1 |«DMX( 2 I 
OMXI  21 aOMXt  31 
0MXl31=rMX(4| 

220  CMX(4)«0 

IF|  KHATUll,GE,XhAT(  JZMCa  TC  226 
IRLMXaO 

225  |F( XHAT < J4t •GF,XHAT ( J3) IGO  TC  227 
IBPMXaO 
227  J1-0MN(1I 
J2-0MN( 2) 

J3bDHN( 31 
J4-0MN(  4 I 
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J 'C  C -trl  1 » = ^ 

IM  XK.-.  T < J'  t .1.  - ,K<A*  < Jt  » »r.T  ’■r  'J5 
I ’lMN-  1 

^ i“(  XM\T ( j 1 ,1.'-,  xhat  ( j 1 > irn  •'r  no 

I >1.  VIN3  1 
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C=(TFM<  J MTPM(’»»Trv<’)  » /r>ENCM 
Wt  * A ( 

0»  =D  >i{2  ) 

1="  ( X(  1 » .GT,X(  ? I IGO  TC  ?a? 

W?  = X(’  > 
r ? =D  X ( J ) 

G ) TO  2 ^3 
2fiZ  ►2rxC’) 

D2*D  X(  ’ ) 

28  3 s-m  )**2-M?*(B-02)**2»/(M2-KlJ 

IF  < TPm.->,lF  . 0»  0)  GD  TC  2'^^ 

ALT=SaaTC  TEMP) 

SM«x'’  = TP|F*M1«(ALT*«2T(P-01)«*2)/(2*0*ALTI 
£MXR=»  JLD  ASMXRF  *NE1»*SMAX3 
265  T'^MP=  ( R-PP  ) ••? 

S5OP=«0LD  *SS0PF  KNEW*! FMP 

IF  I SSOR .GT.AOOTO. 0) SSCR*4  3 0T0.C 

IF  (TFMP.GT.DSOI  GO  TO  2pe 

C LPaAT‘=^  R TGHTC  ST  AR6n  ARD)  PCSITICN  ANC  VELDC  ITT  ESTIMATE  STATES 
TrMP=0-PP 
VQ3VR4.CVATFMP 
P0tpp»r X*  TEMP 

29?  IF(  I CLK  .LE.  SECO  IGO  TP  SCO 
IF!  I CL<  .GT.SECC  IGCt  TO  306 
IF(  StMPD.L.T,2.0)GO  TO  30C 
TFmPzSN  PO-SNRO 
IF  ( TFMP.L  E.0.0  I GO  TO  300 
IFCTFmp.LF  .BMPRSTIGO  to  300 
SPRTP  =l 
PR  = P 
VR«0.0 
S^^RBSTstTFMP 
300  SNPPaSNPP 
GO  TO  306 

C CALULATF  mCASORFD  LFFt(POPT)  PCSITICN 
305  £SQR=»OLD*SSQRAMNEW*l00C0«0 
30e  IF|  IC:»O,»gr.0,ANO.SPRTR,E0,0)GC  TC  31  C 
ICPQ'l 

IF  I SNOO.LT.  ?,0  ) GO  TC  310 
ICRO»? 

IF! ICLK .GT.SECC »GO  TC  307 
IFI SSQR .GP»22500» 0) GO  TO  309 
GO  TO  308 

307  |F(SSOR*GE*iOOOO«OIGO  TO  110 
30a  ICRO*3 

ABSPR^ABS (PR ) 

|F( ABSP R. GT* ABS (FLOAT ( I SCST ( PNLO ) I ) • OR • ABSPR .GT » ABS ( FLO AT ( I SDST (M X 
1 PO»  I) 1 GO  TO  31 0 
ICRO*'' 

GO  TO  310 

30S  IF ( SPRTP.EO.I )SNPHST»0.0 
310  CO  315  I«lt3 

|PX>OMN (I  I 

iPox=inPD( IPX) 

ONI t )al SOSTI IPOX) 

X(  I |zXHAT(IPXI->V« 

315  CONTINOF 

IF(  MPMN.EQ.I  .OP«  IPLMN.FO.I  ICO  TO  380 
L«ON (T) 

IF(X(t ) .GE«O.O.OP*X l2)*GF*O.C.nQ.X( 3l,GE«Ct0)GO  TO  3A5 
TFM( I ) eX( I IADNI 1I*(X|2I-X(3)I 


1 

I 
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BEST'AVAIMBIE  copy 


T w(  r ) = x(  ? ) *i:.N( 

T \«(  :’)  = X(  "<)  *Dr4(  1 ) *(  X <l  l-X(?  ) ) 
f'<-K  JM=?  ,c  *(  T 'M(  ’ ) M(  a » ♦TFV(3  ) I 
ip( pt  >00  Tr  "jeir 

r.  I ^ a 0 r = • . ' 

TFM  ( 1 1=  Tr  ><(  I ) a iNC  I > 

’4  3 C-’NTrvtJF 

L = ( T • TO  >4  < ■’)>  /0‘^KFM 

M : = X « • ) 

0’  =oN<,-e  I 

1“^  « X(  » ) .OT,X(  •«  ) ICO  Tr  ^4S 
►'''=X  ( ■ > 
n;aPN( ’ > 

GU  T -1  T 4 > 

■>4  ■■  : =X  ( ’ ) 

r'c=0'J( ) 

t T"^  MP=  ( Ml*  (L  -nj  I ♦*2-V5*{  L-OZ  ) **?  I/IMS-Ml  ) 

!r(  T=‘M'^,l  )r,n  TD  343 

4L  T=  003  T<  TFVO) 

OMIM’=TPir**l*(  OLT**r*(L“01 )**?)/(?»0**LT) 

4MNL  = W:3LC  *SMM.«- 
343  T^'Mf>  = ( l-PL)  **2 

S'=0L  = wnLD*35  0L*«>N'-'W*TirMP 

IF(  SSOL  .G  T,*3!3  3r3.  0)  «3CL=4  0 0CC.O 

!>■  ( T="r>.GT,03QI  on  •'C  735 

C vJ°C*Tc  L''FT(pnuT>  OCSITION  AMO  VFLOCITV  FST  IMATE  STATES 
■’r' 3 T-»<r>=L-PL 

VL=VI.*CV*T“MP 
PL  = 3L»C  X*  Tt  MP 

•f  •=  !>=  ( ICLK  .LC  . 5=03  >00  TO  770 
I'^  ( I Cl<  .GT.sr  CC  >00  ”0  3«A 
I’-T  SNI  l.GT, -2,''  >00  iC  7T0 
Tf-  MP  ='5NLC-''-NLP 
IP(T-mP,G3.'3.C  I G'i  ■'ll  TTO 
1F< T' MP.Gf  . SNLOST >G0  'O  TTO 
epp  TL  = 1 
PL  = L 
VI  =9*' 

SMLHS'=TF  MP 
37  3 £M|  P=3MCn 
G"  T (1  3 n? 

7“C  '‘■0L  = w0Lr>*330L*  «(N=A*1  0900»  0 

IF(  13^1  KE  .N‘-'.9  ) S30C'l»Cil.P*S50P*WN5W*1  OCCC.C 
3h5  IF  ( I Cun.N'-,  0«  AMO.SPPTL.FC.C  IGC  TC  OPO 
icltf; 

|F(SNLa.GT»"?,''  )G0  TP  ■'OC 
(TlOs'' 

IF  < I a<  .GT.GECr  >00  TP  737 
JF«  SOOL.r.F,  ?2500.0>  GP  TP  730 
".F  T T 7 33 

7BT  IFtSOOL.GF, 10099, OOP  TO  ■»  0 0 
3(33  I-LO*’ 

AT3PL*AB3(PL ) 

ir<  A30TL, GT, AGO (FLOAT! rOCSTCMALC) > ) • PP , ABSPL ,0 T, ABS ( FLUA T( |SOST(MX 
1 «■  0 ) > ) ) GP  ' P 3 7 0 
I 'L  0(^.9 
GP  T ' •»  00 

3“?  IF( SPT'L.FO.I )SNL8ST»C.O 
•09  IF<  I CLX  ,LT,SFC‘>0)Gr  TO  *00 

C/LU  fq  I ( XH  AT  .••IS,  ir  N . lOPP.  KP  ) 

*00  S'TUON 
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FAILURE  DETECTION-AND-IDENTIFICATION  SUBROUTINE 


A.  Purposes 

To  monitor  continuously  each  sensor's  performance,  and  auto- 
matically to  identify  and  to  reject  failed  or  poorly  performing  sensors. 

B.  Calling  Sequence 

Call  FDI (XHAT,NS, ISN, lORD.NF) 

C.  Inputs 

Through  variables  in  calling  sequence 

XHAT^  Scaled  sensor  outputs 

NS  Total  number  of  sensors  at  a baseline 

ISN.  Indicates  failed  sensors  at  a baseline 

1 

lORD^  Sensor  locations  at  a baseline  with  respect  to 

the  centerline  of  runway 

NF  Total  number  of  failed  sensors  at  a baseline 

The  remaining  inputs  are  passed  through  labeled  common  blocks 
/INOU/  Specifies  input,  output,  and  scratch  device  numbers 

/INPT3/  Contains  constants  to  filter  mean  and  mean  square 

D.  Outputs 

Through  labeled  common  blocks 

/VALDO/  Contains  filtered  mean  and  filtered  mean  squared 

/FDIO/  Contains  filtered  variance,  total  filtered  mean, 

and  total  filtered  variance. 

E . Subroutine  Length 
26 3g  locations 
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Warning  or  Error  Mesaages 


a) 

"SENSOR  1 

"x"  ON 

(OB, IB) 

EXCEEDS 

MEAN" 

b) 

"SENSOR  « 

"x"  ON 

(OB, IB) 

EXCEEDS 

VARIANCE " 
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k 


if" 

I 

FAILURE  DETECTION-AND-IDENTIFICATION  SUBROUTINE 


Initialize  total  filtered 
mean  and  total  filtered 
variance  to  zero. 


Calculate  filtered  mean, 
filtered  mean  square  and 
filtered  variance  for 
each  senior 

Calculate  total  filtered 
mean  and  total  filtered 
variance 


Find  largest  absolute 
difference  between  total 
filtered  mean  and  filteredl 
mean,  and  store  in  BIGM 
along  with  sensor  number 
in  IBM 

i-i:.  :: : 

Find  largest  absolute 
difference  between  total 
filtered  variance  and 
filtered  variance, and 
store  in  DIGV  along  with 
: sensor  number  in  IBV 

^ 1 . 


-7frrB~IGV>25^ 


_ 1> - . - - 

Record  which  sensor  failed 

1 

. . _ . 

Record  which  sensor  failed  to  do 

to  do  large  mean  and  elimi- 
nate from  system 

large  variance  and  eliminate 
from  ^stem 

:TI: 


RETURN 
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BESI'AVAIWBIE  COPY 


FAILURE  DETECTION-AND-IDENTIFICATION  SUBROUTINE 


SCnwrCTlNE  FOMXMAT.NS.ISKfICBD.NF) 

C-AL  MALL 

DIMrNSIQN  I SN  ( 30  | . ISDSTO-  »,  lORDOO).  METSCK)*  XHATOO).  Z<3C) 
COMMON/INOU/ICP. IPfi . IPU.  lOT. ISl  • IS2  t ITV. IS3> ISA. ISS. IS6 
C0Mvl0N/lNPT3/C0LC.CNeii 
CCMWCN/V  ALOO/XS  130  l.rS(3''  ) 

CCMMON/FOIi?/ZS<3.*  > .mall  .VALL 
MALL=  . 

VALL=  ■.  ; 

EIGM=-I,C25 
UIGV=-l ,=2S 
cn  1 > 1*1 .OS 

I0X= lOPO ( I » 

1F< ISN(IDX).EQ.l >GO  TO  lOO 
Tr MiCNe  **XHAT<  I ) 

TcMS=TEM*XHAT< I > 

3C  xs< n =COLC*XS(  1 )*TEM 
Y3(  I >=CCLO*VSn  )ATEMS 
Z3(  1 l*VSI I I-XSI  I )**2 
FALL  = MALLAXS  < I I 
VALL*VALLAZS( I ) 

I'P  CONTINUE 
NGS=NS-NF 
I*ALL  = MALL/NGS 
VALL*VALL/NGS 
oa  2^  '■  I *I  .NS 
iox=iCfio<  I » 

I<^(  ISNIICXI.EQ.l  )G0  to  20f 
ABGMsASSIMALL-XSI  I I I 
IF(ARGM.lE.HIGM)GO  to  iso 

ISMS  I 

e igmsargm 

IS',  APGV*ABS<  VALL-ZS(  I n 

IF(  ARGV.LE.B  IGVIGC  TO  24..'! 

IHV=I 

eiGVsARGV 

2 4!:  continue 

IF(  niGM.LE.S.''  IGO  TO  25f 

NCsNFAI 

lOXs lORDI IBM  t 

lENI ICX )*1 

MALL=  C FLOAT «NGS  >*MALL-XS(  IBM)  )/<NS-NF ) 

NGS*NS-NF 
XS( IBM»sO.? 

•RI  TEI  ITr.9j4''  » IBM 
25..  IF(BIGV.Le.25.0  AGP  TO  3f  0 
NFsNFAI 
I0X*I0R3(ieV> 

ISN<  IPX  1*1 

vall*<floatings)*vall-zs<  iev»  )/(ns-nf  i 

Nr.S*NS-NF 
ZS( 1BV)S>.0 
XRITE ( 1 TV.9LCI i IBV 
3 'I  CETURN 

9'':0  FORMATUOH  SENSOR  • t2.19h  ON  OB  EXCEEDS  MEAN) 

9'’>''1  FORMAT!  1)H  SENSOR  * I2.23H  ON  OB  EXCEEDS  VARIANCE) 

END 


C 
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CK  SUBROUTINE 


A.  Purpose 

To  determine  the  location  of  the  first  non-failed  sensor  in  the 
sensor  line. 

B . Callincj  Sequence 

Call  CK(IS,ISN,IORD,NS,IANS) 

C.  Inputs 

Through  variables  in  calling  sequence 
IS  Sensor  number 

ISN^  Indicates  failed  sensors  at  a baseline 

lORD^  Sensor  locations  at  a baseline  with  respect  to  the 

centerline  of  runway 

NS  Total  number  of  sensors  at  a baseline 

The  remaining  inputs  are  passed  through  a labeled  common  block 

/INOU/  Specifies  input,  output  and  scratch  device  numbers 

D . Outputs 

Through  variable  in  calling  sequence 

IANS  Sensor  number  of  the  first  non-failed  sensor 

E.  Subroutine  Length 

50g  locations 

F.  Warning  or  Error  Messages 

"ERROR  IN  CK  SUB CHECK  FAILED  SENSORS " 
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SUBROUTINE  TO  LOCATE  FIRST  NON-FAIL 
SENSOR  IN  THE  SENSOR  LINE 


';u‘’'^ruT  iKc  rKds.isN.icaCtNS,  i#ns  » 

''CVVc^/I^CU/I'■^,fo'>,lPt^ICT,  1*1  ,IS2  .T'rY,I‘:3,IS4,IS5#ISf. 
ClMPNrirN  T*NM0  »,  »nBD«  3C1 
nr  ion  J.i*-  ,N< 

IDi<»Tn'=n(j) 

ini^Mirxi.KF.nor  -*0  120 

ICC  rONTtNUP 

wt  T’’e(ny, 10001 

I'-cc  Ff*'**!' -iCPFor  CK  «;LP....r»-ErK  fhiec  5ebscp<....p) 

r F'^USN 

120 

pcTijnN 

ENC 


i 


! 


B-22 


APPENDIX  C 


KALMAN  FILTER-GAIN  EQUATIONS 
FOR  A DOUBLE  INTEGRAL  PLANT 


Consider  the  double  integral  plant  system  described  by  the 
following  state  equations 

Xj  = X2  , (C-1) 

*2  = u , (C-2) 

where  u is  a white  Gaussian  process  noise  with  the  following  statistics 

E[u(t)]  = 0 , E tu(t)u(t+T) ] = q6(T)  , (C-3) 

and  q is  the  strength  of  the  process  noise.  Let  m(t)  be  the  system 
output  given  as 

m(t)  = Xj^(t)  + v(t)  , (C-4) 

where  v(t)  is  a white  Gaussian  measurement  noise  with  statistics 

E[v(t)]  = 0 , E Iv(t) V (t+T) ] = r6(T)  , (C-5) 

with  r defined  as  the  strength  of  the  measurement  noise. 

The  object  of  this  development  is  to  obtain  the  estimator  which 
produces  the  minimum  mean-square  error  in  estimation  of  x^  and 
from  the  measurements  m.  This  end  can  be  readily  attained  by  defining 
the  state  vector 


Xj^  (t) 
X2  (t) 


x(t) 


(C-6) 


I 


and  writing  the  vector  matrix  system  equations  as 


x(t)  = Px(t)  + Gu(t) , 


m(t)  = Hx(t)  + v(t)  , 


where  the  matrices  F,  G and  H are 


0 1 


0 0 


(C-10) 


H = II  0) 


(C-11) 


The  minimum-error  variance  estimator  is  then  given  by  the  Kalman  filter 
(13) 

equations 


x(t)  = Fx(t)  +K(t)(m(t)  - Hx(t)l, 


(C-12) 


P(t)H' 


(C-13) 


where  x(t)  is  the  estimate  of  the  state  vector  x(t) , as  inferred  from 
the  available  measurement  process  m(t),  K(t)  is  the  matrix  of 
Kalman  filter  gains  and  P(t)  is  the  error  covariance  matrix  which 
satisfies  the  Riccati  equation 


p(tl  - FP(t)  * P(t)r^  ^ 


{C-14) 


Now  the  vector  matrix  differential  equation  for  the  estimator 
can  be  expanded  to  obtain  scalar  differential  equations  for  the  estimates 
of  the  individual  elements  of  the  state.  These  are 


Xj^{t)  » X2(t)  + )Cj^(t)Im(t)  - Xj^(t)), 


(C-15) 


Xjit)  “ )c2(t)Im(t)  - Xj(t)l, 


(C-16) 


J 


) 


where  the  filter  gains  kj^(t)  and  k2(t)  are  the  first  and  second  elements 
of  the  gain  matrix  given  in  (C-13) , 


k^(t) 


Pll(t) 


(C-17) 


k2(t)  = 


Pl2(t) 


(C-18) 


with  and  p^2  defined  as  the  elements  of  the  covariance  matrix  P(t). 

Pll(t) 


P(t)  = 


Pl2(t)  P22^<^) 


(C-19) 


From  Eg.  (C-14)  the  differential  equations  for  the  scalar  elements  of 
P(t)  can  be  obtained  as 


Pll(t) 

, Pll 

^Pl2  ■ r 

(C-20) 

Pl2 

_ _ P11P12 

P22  r 

(C-21) 

P22 

2 

= q - ^ 

4 r 

(C-22) 

If  the  filter  is  allowed  to  reach  steady -state  stationary  operation, 

the  covariances  approach  constant  values.  These  can  be  obtained  by 

setting  the  left-hand  sides  of  (C-20,21,22)  to  zero.  Eq.  (C-22)  then  yields 


12 


= (qr)^  , 


(C-23) 


and  applying  this  result  to  Eq.  (C-20)  obtains 


11 


113 

J J 

2 IT  • 


(C-24) 


Thus,  from  Eq.  (C-17)  and  Eq.  (C-18)  the  steady-state  filter  gains  are 


h • 2^  (f)*  . 


(C-25) 


C-3 


(C-26) 


k 


and  hence  the  steady-state  filter  gains  are  completely  determined  by 
the  ratio  of  process-noise  strength  to  measurement -noise  strength. 

At  this  point  it  is  useful  to  return  to  the  differential  equations 
for  the  estimator.  If  stationary  operation  prevails,  then  kj^  and  k2 
are  constant.  Taking  the  Laplace  Transform  of  both  sides  of  Eq.  (C-15) 
and  Eq.  (C-16)  then  yields 

sxj^  = X2  + kj^  [m  - Xj^],  (C-27) 

SX2  = k2[m  - Xj^)  . (C-28) 

X2  can  then  be  eliminated  and  the  input/output  transfer  function 
relating  m to  x^  obtained  as 


m 


kj^s  + k2 


+ kj^s  + k2 


(C-29) 


The  denominator  of  Eq.  (C~29)  indicates  that  the  estimator  is  a second 
order  system  with  undamped  natural  frequency 


(C-30) 


and  damping  ratio 


(C-31) 


Upon  substitution  of  Eq.  (C-25)  and  Eq.  (C-26)  these  quantities  may 
then  be  written  as 


1 


(C-32) 


C-4 


and 


C = -ij-  = 0.707  . (C-33) 

2^ 

Thus, the  estimator  always  has  a damping  ratio  of  .707  and  its  band  pass 
is  the  fourth  root  of  the  ratio  of  process  noise  to  measurement  noise. 

The  estimator  actually  employed  in  vortex  estimation  is  implemented 
in  a digital  computer.  Although  the  filter  equations  developed  here 
are  continuous  differential  equations,  the  results  obtained  are 
applicable  to  the  digital  filter  because  the  sample  rate  employed  in 
the  computer  is  approximately  50  times  faster  than  the  filter  response 
time . 
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APPENDIX  D 


REPORT  OF  INVENTIONS 


Under  this  contract  an  improved  analytic  model  has  been  formulated 
for  predicting  wake  vortex  transport  in  the  terminal  area  environment. 
While  the  objectives  of  the  contract  have  been  successfully  fulfilled, 
no  innovation,  discovery,  or  invention  has  been  made. 
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